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The effort to discover a quantum theory of gravity is motivated by the need to reconcile the 
incompatibility between quantum theory and general relativity. Here, we present an alternative approach 
by constructing a consistent theory of classical gravity coupled to quantum field theory. The dynamics is 
linear in the density matrix, completely positive, and trace preserving, and reduces to Einstein’s theory of 
general relativity in the classical limit. Consequently, the dynamics does not suffer from the pathologies of 
the semiclassical theory based on expectation values. The assumption that general relativity is classical 
necessarily modifies the dynamical laws of quantum mechanics; the theory must be fundamentally 
stochastic in both the metric degrees of freedom and in the quantum matter fields. This breakdown in 
predictability allows it to evade several no-go theorems purporting to forbid classical quantum interactions. 
The measurement postulate of quantum mechanics is not needed; the interaction of the quantum degrees of 
freedom with classical space-time necessarily causes decoherence in the quantum system. We first derive 
the general form of classical quantum dynamics and consider realizations which have as its limit 
deterministic classical Hamiltonian evolution. The formalism is then applied to quantum field theory 
interacting with the classical space-time metric. One can view the classical quantum theory as fundamental 
or as an effective theory useful for computing the backreaction of quantum fields on geometry. We discuss a 


number of open questions from the perspective of both viewpoints. 
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I. MUST WE QUANTIZE GRAVITY? 


The two pillars of modern theoretical physics are general 
relativity, which holds that gravity is the bending of space- 
time by matter, and quantum field theory, which describes 
the matter that lives in that space-time. Yet, they are 
fundamentally inconsistent. Einstein’s equation for gravity 


82G 


GY = 5 Tw (1) 


has a left-hand side which encodes the space-time degrees 
of freedom via the Einstein tensor and is treated classically, 
while on the right-hand side sits the energy-momentum 
tensor encoding the matter degrees of freedom which must 
become an operator 7”” according to quantum theory (we 
henceforth denote quantum operators by placing a on 
them). The widespread belief is that this inconsistency 
should be remedied by finding a quantum theory of gravity 
so that the left-hand side of Eq. (1) is also an operator. Yet, 
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although we have candidates such as string theory, which 
is in its mid-50’s [1], and loop quantum gravity turning just 
over 40 [2-4], a convincing theory of quantum gravity 
remains elusive. 

Another possibility is that gravity should still be treated 
classically, but that Einstein’s equation should be modified. 
There is some sense to this. Space-time, though dynamical, 
can be understood as describing the background on which 
quantum fields of matter interact. Gravity, through the 
equivalence principle, is unique in this regard—no other 
fields can be described as a universal geometry in which 
quantum fields live. Having this classical background 
structure appears crucial to our current understanding of 
quantum mechanics. Quantum field theory allows us to 
predict the result of future measurements on a spacelike 
surface, based on past initial data. The identification of a 
Cauchy surface on which to give these initial data, not to 
mention the identification of a one parameter family of 
future hypersurfaces, is natural when the background 
structure is classical, while it is unclear whether it is possible 
to quantize this causal structure in a way which is indepen- 
dent of these background choices. While this difficulty may 
be merely technical, background-independent approaches 
such as loop quantum gravity, string field theory [5], and 
others [6,7] have yet to demonstrate that this is possible. 
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And while the ultraviolet catastrophe required the quanti- 
zation of electromagnetism, the analogous issue in gravity is 
less well understood [8] and muted by the fact that 
gravitational systems in equilibrium occur in exotic sit- 
uations [9], or would emit gravitational radiation over very 
long timescales [10]. 

Nonetheless, there is no strong case that gravity is 
exceptional, and the overwhelming consensus has been 
that we must quantize it along with all the other fields. In 
part, this is due to a number of no-go theorems and 
arguments over the years [11-21] purporting to require 
the quantization of the gravitational field. A recent no-go 
theorem [21] holds not only for coupling classical mechan- 
ics with quantum theory, but also with modifications to 
quantum mechanics, so-called postquantum theories which 
alter the quantum state space or dynamics [22]. But given 
the difficulties faced in constructing such a theory, it seems 
important to revisit this consensus. The question of why 
space-time might be fundamentally classical in contrast to 
other forces will be considered elsewhere [24], while here 
we merely consider the necessity of quantizing gravity. 
Indeed, not only do we find that it is logically possible to 
have quantum fields interacting with classical space-time, 
we construct such a theory. The no-go theorems on the need 
to quantize gravity implicitly or explicitly assume that the 
coupling of quantum fields to space-time is not stochastic. 
This loophole can be exploited to construct a theory of 
classical space-time which can be acted upon by quantum 
fields. 

Although the construction of a consistent fundamental 
theory is the principal motivation of the present work, two 
other motivations remain even if one believes that the 
fundamental theory of gravity is quantum. The first is to 
better understand the backreaction that quantum fields have 
on space-time. In many cases, such as in the evaporation of 
large black holes, or in inflationary cosmology, we are 
interested in the limit in which space-time can be treated 
classically, while the matter itself, such as the black-hole 
radiation, or vacuum fluctuations, must be treated quan- 
tumly. We would like to compute the backreaction that 
quantum matter produces on space-time, yet we have no 
consistent theory which enables us to do this. The current 
approach to describe gravity classically is to take the 
expectation value of the right-hand side of Eq. (1), leading 
to the semiclassical Einstein equation [25-27] 


Gm = im), @) 
eg 


Although this equation is often used to understand the 
backreaction on space-time, it is valid only when fluctua- 
tions are small, while to understand quantum effects, we are 
precisely interested in the case where quantum fluctuations 
are significant. Taking the semiclassical Einstein equation 
seriously leads to pathological behavior. If a heavy object is 
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in a Statistical mixture of being in two locations (as depicted 
in Fig. 1), a test mass will fall toward the point between the 
two locations, and the object itself will be attracted to the 
place where it might have been. This has been ruled out by 
Page and Geilker through the performance of a gedanken 
experiment in which a lead mass was placed in a position 
that depended on the reading of a Geiger counter situated 
near a radioactive source [28]. It will come as no surprise to 
the reader, that they verified that the dynamics depicted in 
the right panel of Fig. 1 is what occurs. 

This problematic behavior of the semiclassical Einstein 
equation has generally been viewed as a reason why we 
must quantize gravity [31-34]. However, one sees the same 
apparent paradox in a fully quantum or fully classical 
theory of two systems. Take for example the interaction 
Hamiltonian Hyg = Hy ® Hg, and the limit where the M 
and G system are close to a product state with density 
matrix pyc ¥ Gy ® Gc. Then, for as long as the two 
systems remain approximately uncorrelated, Heisenberg’s 
equation of motion for the M system, after tracing out the G 
system gives 


where we use units where A = c = 16a2G = 1. The state of 
the system M also appears to be coupled to an expectation 
value, even though we know that this is not what happens. 
By tracing out the G system, we have lost view of the fact 
that the two systems will get correlated. Once the systems 
become correlated, the reduced dynamics is no longer even 
linear in the reduced state 6 since t'gHghyc # (Hg)éy. 
While this in no way suggests that the semiclassical 
Einstein equation is correct even on average, its pathologies 
have led to a widely held objection to theories which treat 
gravity classically, while it is clear from the above example 
that the pathologies of the semiclassical Einstein equations 
have nothing to do with treating one system classically, and 
everything to do with taking expectation values, and thus 
losing sight of correlations. The same apparent pathologies 
occur if both systems are classical, but we allow for 
probabilistic mixtures of states. The theory presented here 
does not suffer from these issues; the equations of motion 
are linear in the density matrix, just like quantum mechan- 
ics and the classical Liouville equation. In order to under- 
stand the interaction of two systems, whether they are both 
quantum, both classical, or are classical quantum, we need 
to account for the correlations between them, and this is 
more easily done by considering the evolution of a joint 
density matrix or probability density. 

This leads us to the second motivation for the present 
work: Even if one does not regard the theory as funda- 
mental, it provides a sandbox, or toy model, for under- 
standing quantum gravity. Indeed, a natural first step to 
quantizing gravity is to begin with the Liouville equation 
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FIG. 1. The prediction of the semiclassical Einstein equation is depicted in the right panel, compared with what we expect to occur 
based on the statistical interpretation of the density matrix in the left panel. If a planet is in a statistical mixture of being in two possible 
locations (“L” and “R’’), then we expect the gravitational field to become correlated with it. A test particle therefore falls toward one of 
the planets. If one treats the semiclassical Einstein equations as fundamental, the test particle falls down the middle, which is indeed the 
average trajectory of its path in the left-hand panel. The planets are also attracted toward the place they might have been. The systems in 
the linear theory considered here exhibit the behavior depicted in the left-hand panel. The case where we attempt to put the planet in a 
superposition of the two locations is discussed further in Refs. [29,30]. If the gravitational field is suitably different in the two cases, in 
the sense of how distinguishable the field configurations are when noise is added to them, then there will be a limit on how coherent the 
superposition can be [30]. In the case where there is coherent superposition, the initial trajectory of the test mass will be insufficiently 
correlated with the planet’s position to enable a determination of where the planet is. However, after the decoherence time, it may be 
possible to determine the planet’s location from the trajectory of the test mass, since its trajectory will become correlated to the planet’s 
location [29]. This is also what one expects from a fully quantum theory, with the stochasticity of the test mass being due to vacuum 


fluctuations of the metric. 


for general relativity and a probability density p(t) over its 
phase space. After all, the probability density and the 
Liouville equation have a lot in common with the quantum 
density matrix and the Heisenberg equations of motion. As 
we see with the semiclassical Einstein equation, some of 
the conceptual difficulties faced in understanding quantum 
gravity are also found in classical probabilistic theories and 
are resolved by considering instead the Liouville equation. 
The classical limit of the theory we present here should thus 
be close to the Liouville equation for classical Einstein 
gravity introduced as Eq. (103). Many issues which arise in 
quantum gravity, especially in the canonical quantization 
framework, will turn out to have a similar but more 
tractable form here. In particular, the difficulty encountered 
in understanding the role of diffeomorphism invariance 
when one does not have a single 4-geometry or single time 
slicing, as well as the technical difficulty encountered in 
calculating the constraint algebra, also appear in the theory 
considered here. And even in some cases, in the fully 
classical limit of the theory. 

Let us now turn to another widely held belief on why a 
classical theory of gravity is incompatible with super- 
positions and the uncertainty relation. This argument is 
usually attributed to Feynman [12,35] and Aharonov and 
Rohrlich [36] in relation to the gravitational field (see also 
Refs. [14,37—40]). As depicted in Fig. 2, if the gravitational 
field produced by a particle is classical, the field can in 


principle be measured to sufficiently high precision to 
determine the particle’s position without disturbance. This 
would reveal which slit it went through, and would either 
prevent particles from being in superposition or would 
violate the uncertainty principle if the interference experi- 
ment of Fig. 2 is carried out. However, if the classical field 
responds to the quantum system stochastically [43], then 
measuring the classical degrees of freedom will not 
necessarily determine the particle’s quantum state to 
arbitrary precision. This could be because the gravitational 
field reacts in an indeterministic way to the presence of 
matter, or because there is only a probability that the 
gravitational field reacts to the position of the particle 
during the time of the interference experiment [45]. In 
either case, the gravitational degrees of freedom contain 
only partial information about the location of the particle. 
We thus see that previous arguments for requiring the 
quantization of the space-time metric implicitly assume that 
theory is deterministic, and are not a barrier to the theory 
considered here. 

The fact that coupling classical gravity to quantum theory 
necessarily requires stochasticity [49] is particularly com- 
pelling in light of the black-hole information problem 
[51-53], and its sharper version, the Almheiri-Marolf- 
Polchinski-Sully (AMPS) “paradox” [54,55]. Since deter- 
ministic theories appear to require a breakdown of gravity at 
low energies [56], serious consideration to nondeterministic 
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FIG. 2. Consider a variation of the gedanken experiment proposed by Feynman [12,35] and Aharonov [36] in which a massive particle 
travels through two slits as a plane wave. Quantum theory predicts an interference pattern on the screen behind the slits which enables a 
determination of the particle’s momentum. The local Newtonian potential it creates can be measured everywhere. At late times, the state 
of the electromagnetic field far from the slits will be slightly different depending on which slit the particle goes through, but these two 
states need not be orthogonal. The quantum nature of the electromagnetic field prevents an experimenter from distinguishing these two 
nonorthogonal states and while the electromagnetic interaction might lead to some decoherence of the particle, if the interaction is small 
enough the interference pattern remains. On the other hand, if the gravitational field is classical and the interaction deterministic, then 
one could imagine measuring the gravitational field close to the particle and during its flight, without disturbing the field to 
unambiguously determine which path the particle went through. One could therefore never get interference fringes (or the uncertainty 
principle would be violated). Two classical states cannot be different but nonorthogonal. Actually measuring the gravitational field is not 
required since classical systems are undisturbed; even if we cannot actually measure the field to the required accuracy [48], it is 
impossible to write down a pure quantum state of the particle which is correlated with a classical gravitational field determined by the 
particle’s position. Since the theory presented here modifies quantum theory so that the interaction is stochastic, the gravitational field 
need not become unambiguously correlated with the particle’s path and the interference pattern is restored. Probability densities over 


classical states can be different and nonorthogonal. 


theories should be given [59]. However, theories with 
information destruction face imposing obstacles. Banks, 
Peskin, and Susskind (BPS) argued that they lead either to 
violations of locality or anomalous heating [60,61], and 
Coleman argued that they merely result in false decoherence 
corresponding to unknown coupling constants [62]. There 
have been a few attempts to work around the issue raised by 
BPS. Unruh and Wald proposed that if the information 
destruction was at sufficiently high energy, then violations in 
momentum and energy conservation would also be at high 
energy, and hence, not observable in the lab [63]. Models 
which are local, yet violate cluster decomposition [64] were 
proposed in Ref. [65], and ones which have only mild 
violations of energy conservation at the expense of some 
nonlocality were proposed in Ref. [66]. However, none of 
these attempts have been theories of gravity, and once the 
backreaction to space-time is considered and the constraints 
of general relativity taken into account, the situation is 
different enough to warrant reconsidering the issue. The 
class of theories presented here may allow us to resolve the 
black-hole information paradox in favor of information loss. 
We return to these points in the Discussion section. 

We find that the issue of anomalous heating takes on a 
different form in this theory but continues to raise questions. 
These may be resolved simply by regularizing the theory, but 
it is also possible that we need to introduce new physics at a 


distance scale larger than the Planck scale, as needs to happen 
in spontaneous collapse models [67—70]. Further research is 
needed because the theory we present contains a number of 
unspecified couplings and choices. These can be shown to be 
greatly reduced if one demands that the theory should be 
continuous in phase space [71]. However, this does not 
completely fix the theory, and there remains some additional 
freedom which we expect to be constrained by demanding 
the theory be regularizable, and by experimental constraints. 
Whether there is sufficient freedom in the couplings to satisfy 
these demands, as well as ones coming from symmetry 
principles such as diffeomorphism invariance, remains open. 

Having argued that nothing forbids the coupling of 
quantum systems to classical ones, let us turn to attempts 
to construct such a theory. An early attempt to couple toy 
models for classical gravity to quantum mechanics was 
made in Ref. [72] using the Aleksandrov-Gerasimenko 
(AG) bracket [73,74] 


where A is a quantum Hamiltonian which can depend 
on phase-space points z, and {-,-} are Poisson brackets. 
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The object /(z;t), about which we say more in Sec. II, is 
the classical quantum (CQ) state at time t—a positive 
matrix in Hilbert space, and a functional of the classical 
degrees of freedom z normalized such that 


/ dattp(z;t) = 1. (5) 


However, the dynamics the AG bracket generated leads to 
negative probabilities [72,75]. Unless we are prepared to 
modify the Born rule, or equivalently, our interpretation of 
the density matrix, any dynamics must be trace preserving 
(TP) to conserve probability, and be completely positive 
(CP), for probabilities to remain positive. So, although the 
dynamics of Eq. (4) makes a frequent appearance in attempts 
to couple quantum and classical degrees of freedom [76-78] 
(cf. Ref. [79]), and while it may give insight into some 
regimes, it cannot serve as a fundamental theory. Other 
attempts, such as using the Schrédinger-Newton equation 
[80,81] (cf. Refs. [82,83]), suffer from the problem that the 
equations are nonlinear in the density matrix, and thus will 
lead to superluminal signaling [84—86] and a breakdown of 
the statistical mechanical interpretation of the density matrix 
[87]. And yet, experimentalists working in the field of 
quantum control routinely couple quantum and classical 
degrees of freedom together, usually by invoking the 
mysterious measurement postulate of quantum theory. 
This presents another path to treat gravity classically, where 
one imagines the system of interest is fundamentally quan- 
tum, but that observables are restricted or a measurement is 
performed on them [75,88—90]. However, for gravity, this 
would require first quantizing the theory, which might well 
be an impossible task. Another approach is to apply a 
measurement to the quantum system instead, and then couple 
the output to the classical degrees of freedom such as the 
Newtonian potential [91-94]. Both approaches, at least 
initially, invoke the measurement postulate of quantum 
mechanics, or an ad hoc field which produces a stochastic 
collapse of the wave function [95-100], which then sources 
the gravitational field. 

Here lies the essence of the problem: If nature is 
fundamentally quantum, our notion of classicality hinges 
on the measurement postulate. Thus, the emphasis in this 
work is to invert the logic: Assume that since space-time 
describes causal structure and relationships between the 
matter degrees of freedom, it is a priori and fundamentally 
classical. Then, examine the implications of this. One 
outcome will be that the quantum degrees of freedom 
inherit classicality from space-time. By treating the space- 
time metric as a classical degree of freedom, we obtain for 
free not only a theory of quantum matter and gravity, but 
also a theory which produces the gravitationally induced 
collapse of the wave function conjectured by Karolyhazy 
[101], Didsi [102], and Penrose [103,104]. Collapse arises 
as a simple consequence of treating space-time classically. 
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Such an approach can be found in Ref. [105] where a 
variation of Eq. (2) is studied, but it is argued that the 
interaction of a quantum system with a classical measuring 
device necessarily leads to decoherence and consistency of 
the measurement postulate. Likewise, Tilloy and Didsi 
[92,93] point out that although they introduce a sponta- 
neous collapse model to source the Newtonian potential, 
this can be thought of as formal, and it is the consistency of 
hybrid dynamics which requires decoherence of the quan- 
tum matter. Poulin [66] also note that the classical coupling 
leads to decoherence of the quantum system. 

The theory considered here, however, cannot be reduced 
to a linear spontaneous collapse model, because the dynam- 
ics is not Markovian or even linear in the density matrix 
when the gravitational degrees of freedom are traced out 
(although it is linear and Markovian when we include the 
metric degrees of freedom). It also leads not merely to 
decoherence, but to a genuine objective collapse of the wave 
function. In particular, there are realizations of the discrete 
theory considered here, where stochastic jumps in phase 
space can be finite and correlate to Lindblad operators. As a 
result, changes in the quantum state correspond to changes 
in the classical system. Consequently, the quantum state 
remains pure and uniquely determined by the trajectory of 
the classical system [106]. Purity of the quantum system 
conditional on the classical trajectory has since been shown 
to also hold in the continuous realizations under natural 
conditions [29]. 

Before proposing a theory which couples classical general 
relativity with quantum field theory, we first construct a 
general framework for classical quantum dynamics. In order 
to do so, we derive general linear dynamics which maps the 
set of classical quantum states to itself. More precisely, we 
wish to find the most general linear dynamics which is 
completely positive, trace preserving, and preserves the 
division of degrees of freedom into classical and quantum. 
Two examples of such dynamics have previously been 
proposed. In Refs. [46,107], a master equation with diagonal 
coupling to Lindblad operators £,, was considered of the 
form 


— 5 W(z){LIL,, a t)},, (6) 


where z are discrete classical degrees of freedom, and {-}., is 
the anticommutator. When W%(z|z’) = W%(z)6(z — z’), this 
reduces to the Gorini-Kossakowski-Sudarshan-Lindblad 
(GKSL) or Lindblad equation [108,109] controlled by an 
external degree of freedom z. Equation (6) has been shown 
by Poulin to be the most general master equation in the case 
of discrete classical variables and bounded Lindblad oper- 
ators by embedding the classical system in Hilbert space 
[110]. However, for trajectories which are continuous in the 
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classical degrees of freedom, it will turn out that nondiagonal 
couplings are required. A master equation due to the early 
work of Didsi [47] corresponds to the constant force 


case H(q) = f dx q(x)L(x) or sum of such terms [113], 


Aq) =X; f dx qi(x)Li(x) 


Iss D (H(z). pleit)} - i(q).A(z:0) 
+5 ((A(a).Alst)} — (Ale), Aa) 
+5 / dxdx'{q,(x'),{q)(x), DY (x.x)a(es)} 
+5 | dds! Do,;(x,x) 
x |{A(a), p'(x)}. Ble), {A(@), P(e) H]. 7) 


which adds diffusion and a Lindbladian to the Aleksandrov- 
Gerasimenko bracket of Eq. (4), and is completely positive 
provided the matrix inequality D,Do = ; is satisfied. For 
arbitrary Hamiltonians and dynamics which is continuous in 
phase space, one should turn to the most general continuous 
master equation since derived in [71]. Unravellings for 
arbitrary Hamiltonians has since been taken up in [29]. In 
this work, the classical quantum interaction Hamiltonian 
A1(q) is chosen to be minimally coupled, meaning that it is a 
functional of phase-space fields g;(x) and not on the 
conjugate momenta p;(x), while Hy(z) is any purely 
classical Hamiltonian. The constant force case considered 
in Ref. [47] corresponds to an interaction Hamiltonian 
which is linear in a single classical degree of freedom 
q(x), H(q) = f dxq(x)L(x). The constant force case has 
been studied in the context of measurement theory [111], in 
the classical limit of the collisionless Boltzmann equation 
[112], and in a model for Newtonian gravity [113]. The 
discrete dynamics has appeared in discussions around black- 
hole evaporation [66]. 

In Sec. III, we show that Eq. (14) is the most general form 
of time-local [114] classical quantum dynamics. It includes 
the discreet dynamics of Eq. (6) and the continuous 
Hamiltonian dynamics of Eq. (7). Our derivation uses a 
fully classical quantum formalism and gives necessary 
conditions for the dynamics to be completely positive. 
These conditions are also sufficient in the case of discrete 
dynamics, and in Ref. [71] with Sparaciari, Soda and Weller- 
Davies, we derive the necessary and sufficient conditions 
when the dynamics is continuous in phase space. This 
continuous master equation generalizes Ref. [47] to the 
case of arbitrary interaction Hamiltonians and multiple 
Lindblad operators and classical degrees of freedom. 
Since the continuous master equation requires off-diagonal 
Lindblad couplings, this can be shown to imply that the 
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dynamics of Eq. (6) has discrete, finite-sized jumps in phase 
space [71]. 

With the general form of dynamics in hand, the next step 
is to restrict ourselves to dynamics which is generated by a 
Hamiltonian. This allows us to apply the framework to the 
Hamiltonian formalism of general relativity. This follows 
from considering a perturbative expansion of the more 
general dynamics presented as Eq. (17). This expansion can 
be thought of as a CQ version of the Kramers-Moyal 
expansion [115,116] familiar in classical stochastic dynam- 
ics. We show in Sec. IV that continuous and deterministic 
classical evolution emerges in some limit. Away from this 
limit, the quantum system can retain some coherence, but 
complete positivity requires the classical system to undergo 
diffusion or discrete, finite-sized jumps in phase space. A 
discrete realization of this dynamics makes use of a finite 
and noncommuting directional divergence, Eqs. (73)-(75). 
This dynamics corresponds to a discrete version of 
Liouville dynamics on the classical system and the action 
of a dynamical semigroup on the quantum system, thus 
providing a natural generalization of both classical and 
quantum evolution. A parameter z governs how continuous 
and classical the total system is, with the system’s trajectory 
in the classical phase space going from finite-sized hops to 
continuous deterministic dynamics as t > 0. In this limit, 
the quantum system decoheres and becomes classical. For 
larger t (or f), the decoherence times are longer, while the 
jump distance in the classical phase space increases 
resulting in greater diffusion in the system’s trajectory. 
The dynamics of Eq. (17) also includes the fully continuous 
evolution of Eq. (7). 

In Sec. IV B, we present a simple example of the discrete 
dynamics, namely, a quantum spin system with a classical 
position and momentum interacting with a potential. We 
see that if the particle is initially in a superposition of spin 
states, then it eventually collapses to one of them through 
its interaction with the classical degrees of freedom. In 
essence, the classical potential acts like a Stern-Gerlach 
device, so that after some time, the particle’s classical 
degrees of freedom uniquely determine the value of the 
spin. Unlike standard decoherence, the collapse of the 
quantum state happens suddenly, when the system under- 
goes a classical jump in its momentum. At that point, 
monitoring the momentum unambiguously reveals the 
particle’s spin. The quantum state of the system remains 
pure, conditioned on the state of the classical system. There 
is thus no need for the mysterious measurement postulate of 
quantum theory, a postulate which has meant that the 
interpretation of quantum theory has been open to dispute, 
and whose problematic nature has received renewed 
attention [117]. 

In Sec. IV C, we then extend the framework to include a 
hybrid version of quantum field theory. In Sec. V, we apply 
this hybrid formalism to the case of general relativity. We 
work in the Arnowitt-Deser-Misner (ADM) formalism [118], 
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treating the metric as classical, while the matter living on 
space time is quantum fields. In the classical limit, the 
equations of motion reduce to Einstein equations, in the 
sense of a CQ version of Ehrenenfest’s theorem, that after 
taking the expectation value on the quantum fields, the 
equations of motion are close to the semiclassical version 
of the Liouville equation [Eq. (21)] applied to general 
relativity [Eq. (103)]. There are nonetheless some deviations 
from general relativity in the classical limit, because there 
are additional contributions relating to diffusion and 
decoherence, with some of the experimental bounds on these 
terms since given in Ref. [SO]. 

The equations we derive allow one to calculate the 
backreaction of quantum fields on classical space-time. 
This allows us to go beyond the semiclassical Einstein 
equation, and we expect it to prove useful in attempts to 
understand black-hole evaporation or the effect of vacuum 
fluctuations in cosmology. The interaction of the quantum 
system with classical space-time causes the wave function 
to collapse at a rate determined by the system’s energy- 
momentum tensor. The stochastic nature of the interaction 
allows for the possibility that information is destroyed in 
black holes, although the quantum state conditioned on the 
metric degrees of freedom can remain pure. The theory 
suggests a number of possible experiments which could 
falsify or verify it, and we discuss this in Sec. VI. We now 
begin in Sec. II by reviewing the formalism of classical 
quantum dynamics and then discuss how the evolution law 
we derive relates to both classical stochastic dynamics and 
quantum maps. The form of the master equation, Eq. (14), is 
then derived in Sec. III. Since the initial posting of this 
manuscript in 2018 [119], a number of advancements have 
been made. We continue to highlight some of these develop- 
ments while maintaining the integrity of the original draft. 


Il. CONSISTENT CLASSICAL 
QUANTUM DYNAMICS 


The state of a quantum system is represented in a Hilbert 
space 7/5 by a density matrix—a positive trace one matrix 
&(t) at time ft. A discrete classical system is represented by a 
discrete random variable z, and the system has a probability 
p(z;t) of being at z. In the continuous case, we consider a 
classical system to live in phase space, although we could 
just as well consider any other classical sample space. 
Phase space is a 2n-dimensional manifold M of points 
z =(q,p)', and the state of the classical system can be 
represented by a probability density—a non-negative dis- 
tribution p(z;t) defined over z at time f normalized such 
that [ dzp(z;t) = 1. 

Here, as in Refs. [46,47,66,73,74, 107,111,112], we con- 
sider a hybrid classical quantum system which lives on both 
the Hilbert space and the classical phase space, an illustrative 
example being the classical quantum qubit of Eq. (77). The 
state of a hybrid system f(r) at time f is defined over an 
interval around z by a positive matrix A(z; t) in Hs [120], 
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such that the normalization condition Eq. (5) holds. In 
quantum theory, such a system is sometimes said to be ina cq 
state (or CQ state), and has a density matrix of the form 


Pay = f deplest)l (el ® a(e:0) (8) 


with 6(z;f) a normalized density matrix and |z) an ortho- 
normal basis. Since g and p do not commute in quantum 
theory, if we want to encode phase-space degrees of freedom 
in a quantum system, we could consider two separate 
subsystems, one which encodes qg and one which encodes 
p. Here, we do not consider such an embedding, but rather 
describe the classical quantum state as 


P(z5t) = plz; 1)6(z; t) (9) 


with p(z;t) = trA(z;t) a phase-space density and 
J ttp(z;t)dz <1 when integrated over any interval. 
Conditioned on the classical system being in state z, the 
quantum system is in the state 6(z; f), the normalized state of 
the quantum system. In finite dimension, we can also write 
the CQ state as a matrix whose diagonal entries are phase- 
space densities, as with the CQ qubit of Eq. (77). When we 
consider gravity, each point z is a classical field (the 3- 
metric, and 3-momentum) over R?*. 

We now ask, what is the most general dynamical 
map which take any cq state at time f = 0 to another cq 
state at a later time ¢? This requires us to adapt Kraus’s 
representation theorem [121] to the classical quantum case. 
Table I compares the classical-quantum map to the Kraus 
map for quantum systems and the stochastic matrix for 
classical systems. Table II compares the corresponding 
master equations which we now review. Let us recall that 
for quantum systems, a linear map €(6) which is com- 
pletely positive trace preserving at every time f (we say it is 
infinitesimally divisible [122]) is always generated by the 
GKSL or Lindblad equation [108,109] 


OGty taxes 
a? = — F(A), (0) 
— 1°6(1) | £,A(t)L} -F {lpia A()},| (10) 


with A(t) the Hamiltonian, the set {Z,.} a traceless basis of 
operators orthogonal in the Hilbert-Schmidt norm, and 
A°# (t) are elements of a positive-semidefinite (PSD) matrix, 
which we write as A > 0. Here, and throughout, we use 
Einstein summation notation when there are raised and 
lowered indices. One can always diagonalize A, in which 
case £61; is sometimes called the jump term since one 
could interpret it as giving the rate at which the system 
jumps from being in another quantum state into 6, while 
1{ZLiL,,6}, is sometimes called the no-event term and 
acts to preserve the trace. If the map is non-Markovian 
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TABLET. Summary of Secs. I-III of the finite-time map acting on states for stochastic classical, quantum, and 
classical quantum dynamics. We use Einstein summation notation, so an index which appears both raised and 
lowered should be summed over. In the classical case, the probability density p(z; t) over phase space z is given in 
terms of the earlier probability density p(z;0) and the transition probabilities P(z|z’;t) which define a classical 
channel. In the quantum case, the necessary and sufficient condition for complete positivity is that A(t) as a matrix in 
af is positive semidefinite [which we write as A(t) > 0]. {L,} are a basis of operators. The classical quantum 
positivity and norm conditions combine those of the purely classical and purely quantum dynamics, and generalizes 
the Kraus decomposition theorem to the hybrid case in Eq. (34). 


Map Positivity condition Norm condition 
Classical p(z;t) = f dz'P( |e" ; t)p(z', 0) P(z|z’5t) > 0, Vz, z! fdzP(z|z’,t) =1 
Quantum 6(t) = A%(t)£.,6(0) Ly A(t) = 0 A (t) Like =I 
Classical quantum = A(z; 1) = f dz/A%(z\z’; tL pl; 0)L; A(z|z’5t) > 0, Wz, 2’ fan ale DLL, =I 


TABLET. Summary of Secs. [and II of the master equations governing classical, quantum, and classical quantum 
dynamics, along with the positivity condition in the case of continuous dynamics, or discrete jumps in the classical 
degrees of freedom. In the classical case, continuous dynamics corresponds to the Fokker-Planck equation specified 
by a drift term Dj (z) and a positive-semidefinite diffusion matrix D>(z) = 0 with elements D, ;;(z). One can always 
choose Ly = 1, and the remaining Ly to be traceless [71]. The CQ master equation is a natural generalization of the 
classical rate equation (or Fokker-Planck equation in the continuous case) and the Lindblad equation, and is here 
presented as Eq. (14) (in this table, we use the convention of including the 1°? in the W% terms). In general, the 
coupling constants can have an explicit time dependence which we drop for ease of presentation. 


Classical 
Master equation 0,p(z3t) = f dz’W(z\|z')p(z'; t) — f dz’W(z'|z)p(z; t) 
Positivity (jumps) W(z\z)>0, Vz #2’ 
Positivity (continuous) W(z|z') = —D, i(z')0z6(z — 2’) + Da ij(z')dz 07 5(z — 2’), Dz = 0 


Quantum 
Master equation 0,4 = —i[H, 6] + hPL GL} —}{hPL}L,, 3}, 


Positivity condition h®® a PSD matrix in a, B (h = 0) 


Classical quantum 


L,p(z3 1) Li —4 f dw (2 |<){L} Lag, P(z3 t)} 4 
, W PSD in af for z = z' 


Master equation 0,p(z; t) =, dz Ww (z|z’) Zz 
W” PSD in af Vz 4 2! 


Equation (49) 


Positivity (jumps) 
Positivity (continuous) 
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(or more precisely, not infinitesimally divisible), but has an 
inverse, then it can still be represented by Eq. (10), with the 
positive-semidefinite condition on the 4%4(r) lifted, pro- 
vided the finite-time map is still CP [123]. In finite- 
dimensional Hilbert spaces, CP maps have an inverse, 
but this condition can fail in infinite dimension [124,125], 
in which case, it may not be possible to write the master 
equation in time-local form. 

For a classical stochastic process, on the other hand, an 
initial state at tf = 0, p(z;0) evolves into 


p(z:t) = / dzP(z\z':1)p(2';0) (11) 


with P(z|z’st)>0O and f dzP(z\|z';t) =1. This simply 
ensures that any initial probability density p(z;0) remains 


positive and normalized. This is equivalent to P(z|z’;t) 
being the conditional probability that the state makes a 
transition to z at time ¢ given that it was initially at z’. Fora 
time-local process, the most general dynamics is given by 
the rate equation 


ne : = f adcW(clesole'st) - WEeit)ole) (12) 


where W(z|z’; t) are rates for the system to transition from 
state z’ to state z (this can be thought of as the jump term). 
To conserve probability, we require 


W(z;t) = paewe z3t) (13) 
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being the total rate for the system to transition away from 
state z to one of the states z’. W(z; t)p(z) is analogous to the 
no-event term in the GKSL equation. 

The dynamics can be discontinuous in z, containing 
discrete finite-sized jumps, in which case it is necessary and 
sufficient that W(z\z’;t)>0, Vz#2z' for the dynamics 
to be positive. When W(z|z';t) = —D, j(z')d,6(z — 2’) + 
Dy ii(2')07025(Z —z’), the dynamics is continuous in z, 
and positive if and only if D, ;; is a positive-semidefinite 
matrix in ij corresponding to diffusion. Dj, is called the 
drift term and includes Hamiltonian evolution as a special 
case. W(z;t) = 0 by integration by parts in Eq. (13). The 
continuous dynamics is sometimes called the Fokker- 
Planck equation or forward Kolmogorov equation, while 
the general form of the dynamics is sometimes referred to 
as the differential Chapman-Kolmogorov equation [126]. 

Here we present the most general classical quantum 
dynamical map which is linear and completely positive. In 
Sec. HI, we find the generators of this map, deriving 
the most general continuous-time master equation for 
the dynamics. As in the proof of the generality of the 
GKSL equation, the proof presented here is valid for a 
separable Hilbert space, and countable set of bounded 
Lindblad operators Ly, or alternatively, when the map has 
an inverse at infinitesimal short times. In the case of 
unbounded operators [127], one can cast the equation in 
the form we derive, but it is possible the dynamics are not 
unique. 

Under these assumptions, we find that the most general 
master equation can be put into the form 


= ce oe P(z; Pha 


a” 


+ i dW? (z\2/3t)Lp(2s0)L} 
I nian 
— WH (z5t){L pba Pt) }. (14) 


The condition that the master equation preserve probability 
requires 


vote = faowweey 9 


for each a, f. One can verify this is the case by taking the 
trace of Eq. (15) and integrating over phase space, and then 
requiring the total probability to be conserved in time 


0 x . — 
° ‘| detrp(z;t) = 0. (16) 
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This norm-preserving condition is analogous to the purely 
classical case of Eq. (13) and the pure-Lindbladian case. 

Our derivation is general enough to allow for the 
coupling terms to be time dependent, which would allow 
for the inclusion of various non-Markovian effects, but 
following our general derivation, we drop the time depend- 
ence from 4%(z; rt) and W*(z\z’;t) (one refers to such 
dynamics as being autonomous). For the dynamics of 
Eq. (14) to be infinitely divisible (i.e., CP as a generator 
at every time f), it is necessary that W7(z|z’;t) be a 
positive-semidefinite matrix in a,f Vz#z’, and the 
matrix with elements 4% + W%4(z|z;7) must be PSD. 
These are sufficient conditions for discrete, jumping 
dynamics while the conditions on W(z\z’;t) in the 
continuous case are derived in Ref. [71] and previewed 
here as Eq. (49). For continuous quantum systems, the 
sum over af is replaced by an integral, an example 
being the case where we take Lindblad operators to be 
local field operators written in terms of the coordinate x. 
L, can be any operator, but we can always rewrite the 
above equation in terms of a basis of Lindblad operators 
orthogonal in Hilbert-Schmidt norm, such that the £, 
(distinguished by having a in bold font) are traceless, 
and fy = 1. 

The terms proportional to W(z|z’) and W(z) are 
thus purely classical and account for both stochastic and 
deterministic evolution of the classical degrees of freedom. 
In the case of continuous deterministic dynamics, this term 
is just the Poisson bracket {h°°(z), A(z; t)} with h°(z) a 
classical Hamiltonian. In Eq. (14), we separate out the 
Lindbladian term with couplings 474(z) to be consistent 
with conventions of the classical literature, but it is often 
convenient to include them in W%(z|z’; t). The W79(z|z’) 
can be interpreted as the jump term in both phase space or 
Hilbert space or both since it can induce transitions in both 
the classical quantum system. 474(z) causes transitions in 
Hilbert space only, and W°°(z|z’) transitions only in phase 
space. The anticommutator terms are the collective no- 
event terms. 

In the case that the quantum and classical systems do not 
get correlated in z (for example, when one system acts as a 
large heat bath weakly coupled to the other), the dynamics 
derived here reduces to a semiclassical version of the 
stochastic master Eq. (12) when we trace out the quantum 
system, while if we integrate out z we get an averaged 
version of GKSL dynamics. However, more generally, 
tracing out the quantum or classical system results in 
dynamics that are non-Markovian, since there can be 
correlations between the classical and quantum systems, 
and the reduced dynamics are not even linear on the quantum 
state. When we integrate over z, the presence of correlations 
between the classical system and the quantum one translate 
into different effective coupling strengths W(z) for different 
quantum states [128,129]. 
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In Sec. IV, we introduce a class of theories whose 
dynamics is close to Hamiltonian dynamics. This master 
equation can be expanded as 


op(z;t) 


Po HLA) Az) 


ap e - Tunes . 
+—|E Plat), -s{h pla bat} 
~ EXP (2) Vola nhy + oh.V (ane 
; BI). Tri 
+rh,V-D(2)p(zit) VE +e, (17) 
where V = (0,,0,)’ is the divergence in phase space 


Z=4q, p, and x a matrix of vector fields which, as we 
detail, is generated by a Hamiltonian 


Ale) = Qh, (18) 


over quantum and classical degrees of freedom via the 
symplectic matrix 


0 Tn 
O= : 19 
i 0 | ( 
The couplings 
Xi (z) = QVA®(z) (20) 


are the terms which correspond to Hamiltonian evolution in 
the classical limit. y7’(z) is taken to correspond to a friction 
term and is a vector in the 2n-dimensional phase space 
at each z, with components 15 (z)s while the D(z) 


corresponding to diffusion terms define a matrix both in 
ng 
typically considers D%(z) to depend only on g and the 
diffusion to be in p. Here, we include the pure-Lindbladian 
term proportional to 4%4(z) in the W°?(z) term. 

Dropping terms of higher order in + and taking the 
expectation value on both sides by performing the trace 
over the quantum degrees of freedom gives a semiclassical 
version of the Liouville equation 


af and in phase space with elements Df; (z) at each z. One 


dp(z;t) 


= tr{A(z), p(zt)}} +... (21) 


If we retain the term proportional to c, an appropriate 
choice of y(z) and D(z) corresponds to a hybrid version 
of the Fokker-Planck equation. Truncating the expansion in 
Eq. (21) can result in negative probabilities, but the full 
evolution equation is completely positive and provides 
a natural generalization of both classical Hamiltonian 
dynamics and the GKSL equation. 

We also see that taking the trace over the quantum 
system does not result in the pathological behavior found in 


PHYS. REV. X 13, 041040 (2023) 


the semiclassical Einstein equation (2) as depicted on the 
right side of Fig. 1. Instead, the dynamics of Eq. (21) is 
linear in the total density matrix and similar to the situation 
described in our discussion of Eq. (3) where we have linear 
evolution of two coupled systems which account for the 
correlation. For example, the correlated state of two 
components labeled L and R, 


n 1 r 1 , 
P(z) = 5 PL lz)6r Al 5Pr(z)ér. (22) 
would evolve according to Eq. (21) as 


PL) _  t0Fi( 261. px(c)} +5 (WAZ) be. pele} 


ot 2 
(23) 


with the component p; (z) evolving differently into pp(z) as 
in the left side of Fig. 1. One can also think of 0/dp in 
Eq. (21) acting on each matrix element of A(z; t), while the 
velocities @H(z)/dq are a matrix, and thus the backreaction 
force depends on which state the quantum system is in and 
correlates the quantum system with the classical one. 


I. THE CLASSICAL QUANTUM MAP 
AND MASTER EQUATION 


We first derive the classical quantum map which takes a 
cq state f(z; 0) to another cq state f(z; f) and then derive the 
continuous-time master equation. Recall that we restrict 
ourselves to the case of trace-class operations. We also 
restrict ourselves to the case where the map has an inverse 
at infinitesimal times (which need not be completely 
positive), so that we may write the master equation in 
time-local form [123]. Recall also Kraus’s theorem, that the 
most general map on quantum systems €(6) that is CP and 
TP can be written so that 


(1) = E(6(0)) = )K,6(0)K; (24) 


with the operators K,, satisfying 
KiK, =1, (25) 


where CP means that 1 @ € acting on Hs ® Hg takes 
positive operators to positive operators, and together with 
the trace-preserving condition, ensures that density matri- 
ces are mapped to density matrices, even if the map acts on 
part of a system. 

Let us now consider the case of a map which acts on cq 
states. Because the map has to be linear in A(z; t), we can 
write the total map A, acting on /(0) in terms of quantum 
maps A,)..., such that 


plait) = f dz Azo.Al2'0). (26) 
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Since the initial state could be A(z’;0) = 6(z’ — z9)6(z’), 
Le., a normalized quantum state, with vanishing support a 
finite distance away from some initial point zg, the map 
A.jz. must be completely positive for all z, 2’. The initial 
state which could be peaked around any z’ must be mapped 
to another positive (but not necessarily normalized) quan- 
tum state around z. In the case of discrete z, the maps need 
not preserve the trace; they must merely be completely 
positive, trace nonincreasing, since if we act on the 
normalized state around some point z’, a final state around 
z cannot have a trace larger than 1, as this would correspond 
to a probability f dzp(z;t) larger than 1. Thus, one easily 
sees that a necessary and sufficient condition for a cq map is 
that each A,.,, must be a completely positive, trace 
nonincreasing map, with the total map 


Nyt = f dedye. (27) 


being trace preserving. This is the analog of the classical 
condition { dzP(z|z’) = 1. In the case where z is continu- 
ous, the above statements hold, except that because trA(z; t) 
is a distribution, we should interpret statements about the 
nonincreasing of its trace as referring to its trace smeared 
over an indication function [ dztrp(z; t)é(z). 

We can then write the map in terms of sets of Kraus 
operators K,,(z|z’; 1), 


plait) = f del OK (cles )A(e0)K, (2k (28) 


where 5°, K,(z\z’3t)'K,(z|z/;t) <1. This follows from 
the fact that the action of the entire map must preserve the 
trace; thus, summing over y and taking the trace of both 
sides of Eq. (28) on an arbitrary initial state A(z’;0), 
requires 


Df aekylcles)*Ky(eke'st) =1. 29) 
H 


Equation (28), together with the normalization constraint of 
Eq. (29), gives the most general map on hybrid classical 
quantum states. It is the hybrid generalization of both the 
classical probability map of Eq. (11), and in some sense, an 
extension of the Kraus representation theorem Eq. (24). 
However, it can also be considered a special case of a Kraus 
decomposition, as we could restrict ourselves to quantum 
states that have the form of Eq. (8) and write its action as 


Dealt) = ‘| dedz!S“|z)(z!| ® Ky(<le's1)Peq(O)Ie’) (cl 
@.K,(zlz'32)". (30) 


where the set of Kraus operators is {|z)(z'| @ K,(z|z’;1)}. 
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We now want to derive the generator £ of this map 


A, = exp ( ["ae.) (31) 


n : z3t+6t) —p(zt 
Lplzst) = tim 7 cee. (32) 


using 


A 


This is possible, if at short times, the map has an inverse 
(which need not be CP), since then it can always be 
written in time-local form [123] (the proof contained in 
Refs. [124,125] also applies to maps on CQ states). Le., 


p(t) = K(n)p('s0). (33) 


To this end, we write the Kraus operators in terms of an 
orthonormal basis of operators Ts to obtain 


plait) = > favattelesneaeoe}, (4) 


a,p=0 


where complete positivity for all z, z’ requires the 
A“? (z\z/;t) be a positive-semidefinite Hermitian matrix 
in a, f for all z, z’ (see Ref. [130] and also Appendix A). 
The trace-preserving condition Eq. (29) requires 


y i ded (c\z:1)E1L,, =1. (35) 


ap 


Since {£,} is an arbitrary orthogonal basis of the 
operator space, we can always take £) = 1, and the 
orthogonality condition tra = 0 for a # P means that 


the other L, are traceless, and we denote them with bold 
font a. We are interested in changes of the state over 
infinitesimally short times dt. Since the dynamics is time 
local [Eq. (33)], we can expand the A? as 


NOO(Z|z5 t+ St) = AM(|2/5 1) — H(z — Z)YP(Z's 1) 5t 
+ W(z\z/; t)ét + O61’), (36) 


AP (z|2/3t + St) = AP (z|2'3 2) + 6(z — ZAP (z'; t) dt 
+ W*(z|z’;t)6t + O(St7), a6 #00. 
(37) 


This is just a Taylor expansion where we divide up the terms 
proportional to d¢ into two contributions for later conven- 
ience and for comparison with existing convention. If the 
map acts on an initial state p(z;r), then A(z|z’;t) = 
16(z — z’), since when 6t > 0 the system has to remain at 
the same point in phase space and the quantum system does 
not change. The other A%4(z|z’; t) = 0, since these terms are 
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associated with Lindblad operators Ly acting on the state and 
nothing happens to the state in the limit 6¢ > 0. This leaves 


A(z 2'; t+ St) = 6(z — 2')(1 — dty(z'; 1) 
+ W(z\z'; t)6t + O(dt7), (38) 

= 6(z—2')A%? (z's 1) dt 

+ W*(z\z';t)5t+O(650’), 


A*®(z|z/;t+6t) 
ap#00. (39) 


The first term proportional to df in Eq. (39) can be 
thought of as the probability per time dt the system does not 
move in phase space while the quantum system has 
Leplz tL} applied to it, while the W(z|z’;r) term in 
both equations gives the rate that the system goes from z’ to 
zand has L,p(z; tL applied to it. We adopt the convention 


that the sign in front of 7 is negative, anticipating that it 
will contribute negatively; it is the no-event term, and the 
probability that the system is unchanged can decrease only 
with time. 

The expansion for the A° component is similar to the 
one used to derive Eq. (12) for continuous classical 
stochastic processes. Taking the limit dt > 0 and keeping 
only the first-order terms usually require that the set of 
operators £,, be bounded in order for such terms to be 
small, but here we replace that typical assumption with the 
assumption that the map is invertible at short times [leading 
to Eq. (33)]. To find the generator of the dynamics, we 
substitute Eqs. (38) and (39) into Eq. (34) to get 


plait + 61) = f drole 2)((1 = ayes) AC's) 
+ tA (z; t)Lgp(z's t) + 5tAP (z; t)p(z'; t)L} 
+ 614% (2; 1) p(z'; t)L}) 


si > f acarwer cl tL (x)p(z / t)L 


ap 


(40) 


and subtracting A(z; t), dividing by 6¢, and taking the limit 
dt > 0 gives 


Op (z;t) = —/(z; 


a y(z; t)p(z;t) + 1% (z; t)Lep(z; t) 


+ AB (z; 1) pz; t) Lf + A(z; 1) Ladle; 1) Lf 


+ > fa W (z\2/; t)L,p(z's tL}. (41) 
ap 


As in the derivation of the GKSL equation, we now define 
the two Hermitian operators 


a 1 e i 
A(z) = > ("Lg — aL), 


A 


1 x in 
Mash —yMe) + CAME} + AL), (42) 
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so that we can write Eq. (41) as 


Op(z; t) = ~i[A(z; 


at 1), B(zs)] + 2% (23) A(z Ladle: Lj 


+ f awn clenLpeint, 
+ {N(z51), A(z t)},. 8) 


Having subtracted f(z; t) from Eq. (35) the trace-preser- 
vation condition Eq. (27) becomes Eq. (16), i.e., 


ip (z; t) 
faa Pee ai = 0 


for any input state A(z’; t), which requires 
1 ee 
s)=9 (Soe pla 
op 
+ / dz) -W®(z\z'; HhjL). (44) 
ap 


Defining the outgoing flux W?’(z; 1) by Eq. (15) gives the 
desired form of the master equation (14). 

The first term of Eq. (14) is the free evolution of the 
quantum system at each point in phase space, the second 
term corresponding to 4%°(z; rt) is pure Lindbladian, with a 
rate determined by the classical degrees of freedom. The 
remaining term corresponding to W?? gives the interaction 
term between the classical and quantum degrees of free- 
dom. It is worth noting that the term corresponding to 


ap = 00, 


Vlei) = > dz'W(z\z!)p(z'st) — W%(z)A(zst), (45) 


is the purely classical evolution of the system and is 
identical to Eq. (12). It can be stochastic or it can 
correspond to deterministic dynamics generated by a 
classical Hamiltonian H(z), in which case, we can write 
it in terms of the Poisson bracket, or equivalently, as 


Wz; 4) = XP. Vp(zi0) (46) 
where X°° is the Hamiltonian vector field 


x0 = OVA. (47) 


Recalling Eqs. (20) and (18) 


One _ ah \* 
xp = (5-5) (48) 


V could also include a diffusion term as in the Fokker- 
Planck equation, but it will be natural to treat it in the same 
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way as the other W?, in which case, it will be a stochastic 
version of Hamiltonian dynamics. 

Let us now return to Eqs. (38) and (39) and the condition 
that complete positivity places on A%*(z|z’; t). In general, 
all that is required is that the total map from some initial 
state at t = 0 to some final time is completely positive. In 
this case, the generators at intermediate times t need not 
generate completely positive dynamics if one were to 
replace the current state of the system f(z; t) with some 
other state as long as they generate CP dynamics on the 
state which had been evolving since t = 0. In this way, one 
can consider non-Markovian dynamics. However, in the 
case where we demand that the generator is infinitesimally 
divisible, i.e., generate CP dynamics at all times on any 
state, then the positivity conditions for A?(z|z’; t) to be 
completely positive can be determined from the right-hand 
side of Eqs. (38) and (39). For z#2z’ [we can smear 

| 


n!} 


W(2|z’3t) = 


1 M(t) VO"5(z 


n=0 n! 
with 


M, = M,M;'Mi, (I—M )M>')M, =0. (50) 


Here, Mj! is the generalized inverse of the PSD matrix with 
elements Kg ,M, isa matrix in both a, i where i indexes the 
components of z. It has entries M{%. M9° is a PSD matrix in i, 
j with entries M9 ;- This generalizes the constant force 
dynamics of Ref. [47], and gives Eq. (7), allowing for any 
Hamiltonian A(z; t) to generate classical quantum dynam- 
ics, as well as more general non-Hamiltonian backreaction, 
such as those which include friction. 

In the remainder of the article, we restrict our attention to 
the autonomous case where the coupling constants have no 
explicit time dependence. 


IV. NONCOMMUTING 
HAMILTONIAN DYNAMICS 


We would like stochastic hybrid dynamics, which has 
continuous Hamiltonian evolution in the classical limit. For 
the continuous dynamics of Eq. (49), this is given as 
Eq. (7). In this section, we explore the discrete master 
equation corresponding to Eq. (17) whose first few terms in 
its expansion correspond to classical quantum versions of 
the Liouville equation at first order and Fokker-Planck 
equation at second next order. We begin by taking z to be 
phase-space variables q, p, and later present the case where 
z represents canonically conjugate pairs of field degrees of 
freedom over x, since the generalization is straightforward, 


same M°°(z!,t) VE"5(z 
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Eqs. (38) and (39) over an indicator function with support 
on or away from z = z’ as is done in [71]], A?’ (z|z’;t + 6t) 
is equal to W(z\z’;1)dt, and so it is necessary that 
W*(z|z’;t) be a positive-semidefinite matrix in af. 
When z = z’, we just have the usual GKSL dynamics, 
and so we require W% (z|z’; t) + A%(z, t) be positive semi- 
definite. These condition are also sufficient in the case of 
discrete dynamics, because these two cases completely 
specify the dynamics, and the only condition it has to 
satisfy is complete positivity and norm preserving. The 
necessary and sufficient conditions on W(z|z’;t) for 
continuous evolution are derived in Ref. [71]. They are 
best given in terms of the moment expansion of 
W(z|z’;t), which is not derived until Sec. IV, but for 
completeness we conclude this section with it, expressed as 
a matrix in a, p: 


op 
Zz), Che Me VE"5(z — Z/) 


n!} 


(49) 
z), MP (z,8(z-z’) 


while the notation is more cumbersome. The dynamics is 
generated by a Hamiltonian operator, which can also 
depend on phase-space variables z and a decomposition 
of the Hamiltonian in terms of Lindblad operators L,, given 
by Eq. (18), which we recall as 


A(z) = h® (z)E}E,. 


When we get to the field theory case, we see that Lorentz 
invariance, and later, diffeomorphism invariance, place 
very tight constraints on this decomposition. The h°°(z) 
term is purely classical, since it couples only via 1 to the 
quantum degrees of freedom. There is some ambiguity in 
what is included in the £,, and what in h®*(z), but it will 
turn out that this ambiguity is inconsequential when 
applying this formalism to general relativity. The matrix 
h°? is Hermitian, and it is sometimes convenient to take it to 
be positive as is the case for the coupling term of general 
relativity. We define a generator of dynamics X;, which has 
both phase-space components and components af in 
Hilbert space determined by the Lindblad operators, i.e., 


yet — oh _ on . 

h f,) p ? fs) q 
Just as the Hamiltonian vector field determines the rate 
of flow along different directions in phase space in the 


(51) 


case of purely deterministic classical evolution, here, ew ; 
additionally sets the rate of quantum jumps along different 
Lindblad operators. 

We have already expanded A, in terms of infinitesimally 
small times dt, and we now also expand it in terms of the 
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phase-space vector A :=z-— 2’, which we can then take 
arbitrarily small if we choose. To this end, let us write 
A (z|z/; t) = A®?(z|z — A; t) and perform both an expan- 
sion in 6f and a Kramers-Moyal expansion [115,116] in A. 
We can either do this directly on Eqs. (38) and (39), but it is 
slightly less cumbersome to perform it on W°4(z|z — A; 1), 


Were Asto ge As7) 
= W%(z+A —Alz— A; n)p(z— A; 2) 


7 s ae [We(z+A 


1=0 


atplzt)|. (52) 


We now integrate this expression over A to write 


jaawer(al — A; t)p(z — A;t) 


in terms of the moments 
Mi et = [aatayorwe( + Alz;t). (54) 


Here, we write M?(z; t) as an n-fold tensor power, since 
A®”. V8" = (A-V)" with A-V the directional diver- 
gence in the direction A. This has the advantage that the 
notation is compact, but the disadvantage that it distin- 
guishes the moment corresponding to 0,0, to that of 0,0, 
(for example). The zeroth moments 


Me = / dW (2+ Alzst) (55) 


are just W*(z;t) of Eq. (15). The first moments are the 
vectors 


Me — / dAAW#(z+A 


zt). (56) 


For a Hamiltonian H and decomposition h”’, we see that 
there is a natural choice of W??(z + A|z;1) so that the 
vector of first moments has a contribution corresponding to 
Hamiltonian flow 


X#P(z) = QVh® (2), (57) 


where we indicate that W*’(z + A|z) is generated from the 
field h*’, with the subscript h. This we do in Sec. IVA. We 
can also consider a contribution corresponding to friction 


ye so that 


/ dhbAW(z + Alz) = X%(z) —7*(z). (58) 
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Although there could be some arbitrariness in this division, 

a standard friction term is of the form 7*"(z) = tyah"4/ap 

with y a constant and 7 included for later convenience. 
Next, there is the matrix of diffusion terms 


1 
D(z) = Al dAA ® AW}, (z te Alz). (59) 


Making the Markovian approximation so that we may 
drop the explicit dependence on time, and putting this 
together gives 


jf aawer(ele~ a)a(e- ass) 
= Wy (z)p(zs t) — X9? - Va(zst) + V vp alzst) 
+V-D¥®(z)p(zit)- Vb, (60) 


where we use the fact that V - X Ki = 0, and the notation 7 is 
to indicate that V acts on the second A in the tensor product of 
D(z), while V acts on the first. One typically refers to the 
terms which contain V/(z; 1) as the drift, since they deter- 
mine the rate at which the expectation values of observables 
change. The terms corresponding to second derivatives lead 
to diffusion in phase space. Note that if D?’(z) depends on 
both the g and p, then it can also include contributions to the 
drift, sometimes referred to as anomolous drift. Here we take 
D*(z) todepend only on the g so that there are no anomalous 
drift terms coming from this and higher-order terms. 

The Pawula theorem [131] says that either the Kramers- 
Moyal expansion truncates at the second moment or requires 
an infinite sum of terms. In the latter case, our truncated 
master equation will not be completely positive, and should 
be regarded as an approximation of the dynamics of Eq. (14). 
We call a particular expansion of A in terms of Lindblad 
operators, and the choice of W7?(z + A|z) or A%(z + Alz)a 
realization, and although we present two simple realizations 
in Sec. IVA, the main requirement we demand is that the 
moment expansion include X . (z), the Hamiltonian vector 
field for h’ (z) in its first moment, Eq. (58). This ensures that 
we are able to reproduce the deterministic classical equa- 
tions of motion in the appropriate limit. When we consider 
applying our results to field theory and later general 
relativity, Lorentz covariance and parameter invariance 
implies further constraints on W%(z|z’). 

Substituting the Kramers-Moyal expansion of W7?(z|z’) 
into Eq. (14) gives 


op(z;t) 


5 = HA). Pz] + A%() 


Aa a “A 1 “a “a A 
+ Mq"(2)) |LaPlest) by — 5 {Epa Plet)}, 


— EXP (2) Valzst) bp + La (alone; 


~ : ae 
+£,V-D¥?(z)p(z;t)- Vip 


041040-14 


A POSTQUANTUM THEORY OF CLASSICAL GRAVITY? 


When we consider realizations of W%’, we typically find 
that each order in the expansion carries with it higher 
powers of some constant 7. If we identify (1/7) W%4(z) = 
A(z) + MC? and D%? > cD%, this is Eq. (17) highlighted 
in the Sec. I. 

An important category of realizations are ones in which 
M4? (z) for n > 2 and y%"(z) depend only on the conjugate 
coordinates g, and the diffusion is only in the conjugate 
momenta p, i.e., 


Op(z;t ix 
PSO _ _ A(z), ales) 
WE) [x 1s 
+O Ti alent} —jljte an}. 


~ O ap A 
+ Oey vn (z) A(z t)|L, 
ap(.\> ) rit 
lig Dal apt) ap ett, (61) 


This means that the q still obeys the equation of motion 


0h 
ae (62) 
and allows one to define p deterministically in terms of g 
and g. This can be verified using the Heisenberg repre- 
sentation for all moments of g (see Appendix B). 
Equation (61) is thus a classical quantum version of the 
Wiener process or Fokker-Planck equation, and we refer to 
such a master equation as Brownian. For such a master 
equation, it is the conjugate momentum which is influenced 
by the quantum system. In the more general case, p is 
defined only in terms of expectation values. 

It is worth noting that in our expansion, we have a purely 
Lindbladian term which is itself trace preserving, and a 


total divergence term, with the quantity X*? — y*? — pe . 


Vv acting as a current. Thus, if we integrate over phase 
space, this term also conserves probability, contributing 
only a boundary term. Generally though, this term is not 
positive by itself, and one typically has to include the full 
expansion to ensure this. Looking at each of the terms in 
Eq. (17), we see that the first term is just the deterministic 

| 
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quantum evolution. One can have a purely deterministic 
classical evolution term as well. Since Ly = 1, the X?° term 
is generated by the pure classical term h°° in the expansion 
of A(z) in Eq. (18), 


—X(z) Vr = {h(z), (zi 1)}.- (63) 


The other deterministic term, at least on the classical system 
and for a particular af, is L,X*"(z) - Va(z; HE), which, if 
we take the trace over the quantum system, gives the 
required dynamics we are looking for, Eq. (21). 

The other nondeterministic parts include the pure- 
Lindblad term, and if we trace out the quantum system 
it contributes nothing to the dynamics of the classical 
system, except indirectly, as it evolves the quantum system 
which backreacts on the classical system. On the quantum 
system, it can generate decoherence, and if the Lindblad 
operators are diagonal in raising and lowering operators, it 
raises and lowers the state of the quantum system. The final 
two terms, if we trace out the quantum degrees of freedom 
are the friction and diffusion terms on phase space 
reminiscent of the Fokker-Planck equation, although here, 
there is again some non-Markovianity as the quantum 
system can serve as a memory for the classical system. 

The Lindbladian term can determine the rate at which the 
wave function decoheres, and if this rate is small then the 
diffusion term tends to be large. This is to be expected since 
as discussed in Sec. I, a purely deterministic evolution on 
phase space would decohere the wave function instantly, 
and so a slow rate of decoherence requires either a very 
noisy stochastic evolution or a very slow jump rate. 

It may be worth noting that we can now write a formal 
solution of Eq. (17). One way to do this is to double the 
Hilbert space [132-134] to write the density matrix A(z; t) 
in its eigenbasis 


Alzst) = Dopilz)li) (i (64) 
as a pure entangled state in the Schmidt basis 


A(t) ap = > V Pilz) lida @ |i)e- (65) 


The CQ master equation of Eq. (61) can then be written as 


P i ee ae ee 1 or 
= |-iH(z) @1+i1 @ A(z) + AVY (z) + W%(z)| @ @ L;- 5 hgh @1-51@ it; 


— Ly @ L4XiP Vt Ly @ LBV -1iP(z) + La @ LGV - D¥(z) +--+] |P(z5t)) an. (66) 
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where we use the fact that applying an operator to one 
subsystem of the entangled state is equivalent to applying 
the transpose to the other subsystem [135]. If we designate 
the operator in square brackets as £, then the formal 
solution to the rate equation is 


lA(z5t) ap = e*"|A(z;0)). 


The generator on the doubled system is in some sense a 
Hamiltonian, albeit one with an imaginary part. 


(67) 


A. Realizations of W(z,z—- A) 


Here, we present some simple and constructive examples 
for W%(z,z—A) which have the properties we require, 
namely, positive-semidefinite matrices whose first moment 
contain a term which generates Hamiltonian flow via 
Eq. (58). A simple realization valid for small 7 is 


WP(z,2— A) BA —eXH(e)]. (68) 
This gives the Liouville term as the first moment as 
required, and 1/z as the zeroth moment governing the 
collapse rate. We could take t to depend on a, /, or include 
a distribution f% over z. The realization of Eq. (68) can be 
considered an approximation of 


1 i 
/ dW (z,z— A)P(z— Ast) =— eI A(z; 1) (69) 


with the right-hand side being the exponential of 
the Poisson bracket evolution {h%(q,p),-}:= 
dxl(ah”? /0q,)(0/Ap,) — (0h? /Ap,)(0/0q,)| which gen- 
erates a finite shift in phase space along the flow generated 
by A (dropping the af index for simplicity). This is the 
classical analog of the evolution operator induced by the 
-iftt 


quantum Hamiltonian H, e since 


« Oh 7 oh 
ett-}a(g, p:0 =o( -| di —., [at 0) 
P(g. P30) = Pl 4 me op dn oa 


(70) 


Then e7{""} — 1 acts like a finite-time divergence along the 
direction of Hamiltonian flow generated by h. 

This leads to a natural noncommuting generalization of 
this directional divergence which for a CQ Hamiltonian 
with a single Lindblad operator H =hAL*L, and for 
small z, is 


Wig. p) =—eMILEPL- = {LIL a}, (71) 
In the case of many Lindblad operators, let us define 
basis vectors e% and define 
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1 
Vopb(eite = =| f dAWiP ole A)L Ble Asn} 


=Wr @)5 {Eile Bla t)}4 (72) 


leading to 


Let I p(z thy — 3 {Lpha P(e t) hy 


T 


VupP(z; ne" ~ 


(73) 


which for small ¢ looks very much like a finite-directional 
divergence (both along the direction x7, i 
Lindblad operators). 

This definition depends on the decomposition of Al into 
Lindblad operators, but one can first diagonalize h*/(z), in 
which case, we write the left-hand side of Eq. (73) as 


V4,P(z;t). In cases where 


and along the 


a 1]. A apy rt 1 page A 
V.pp(z; the? x — | Lap(z—tX4? )Lg- 5h pha bls t)},], 


aq 


(74) 


it is sometimes the case that x is non-negative (or 


nonpositive), in which case this discrete: noncommuting 
directional derivative takes the simpler form of a non- 
commuting finite difference 


ap A 
xe VipP (2; t) x 


Care must be taken to ensure it is norm preserving, which is 
the case if, for example, xe depends only on g and the 
jump is in p. Although Eq. (73) looks like a discrete CQ 
Poisson bracket, V,/(z; ft) is only linear in A(z; t) and not 
HZ. On the other hand, the right-hand side of Eq. (4), which 
uses the Aleksandrov-Gerasimenko bracket, is bilinear and 
so it does have some properties more akin to a novel 
bracket although its action is not CP. The same holds true if 
we replace the Aleksandrov-Gerasimenko bracket in 
Eq. (4) with the zeroth and first moment of Vy/(z;t). In 
some cases, we can have a bilinear map which is also CP as 
in Eq. (75). 

The noncommuting divergence is trace preserving and 
completely positive, and allows us to write the master 
equation as 
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If we make the model more deterministic (t > 0), the 
strength of the pure-Lindbladian term increases to com- 
pensate, while if we make the Lindbladian term small, the 
diffusion terms at second order become significant. Here, 
we have temporarily switched units to show fh, since it 
strongly suggests that we set t = h. Thus, as h — 0, the 
decoherence renders the system completely classical, the 
diffusion terms and those of higher order go to zero, and we 
are left with purely continuous classical equations of 
motion. For gravity, the realizations are more tightly 
constrained due to diffeomorphism invariance, but similar 
realizations can be used to obtain dynamics which recover 
general relativity in the classical limit. 


B. Example: Qubit coupled to a classical potential 


Let us now review a simple system solved in Ref. [106] 
with Sparaciari, Soda and Weller-Davies, which contains 
some features we will encounter with gravity [136]. While 
the model bears some similarity to the Stern-Gerlach 
experiment (cf. also Ref. [75]), the behavior of the spin 
can differ depending on how one decomposes the 
Hamiltonian in terms of Lindblad operators. Our system 
consists of a classical particle which couples to the spin of a 
two-level system or the energy of a quantum harmonic 
oscillator. The classical degrees of freedom are the position 


Op(q. P) 


= -i[f(4, p), (a. p)] + {H,Ala. P)} + 


ap 
a 


(0) 


0 


m 


+<[b)(UA P— tBN) HI -{1,2@ P)}y- 


In terms of the matrix elements of the cq state, Eq. (78), the 
master equation decouples into two equations for the diagonal 
components and one for the off-diagonal coherence term 


duy(q.P) pp 9us(4.P) 


= [uy (q.P + TB) —ur(q, p)], 


Ot m oq 
ouy(q.P) _ pouy(g.p)  @ 
=- —7rB)-— 


(79) 


©) .ata.r)] +{Z alan) + 
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q = x and momentum p of the particle, with free evolution 
given by h°° = p?/2. We could take the purely classical 
evolution from W°°(z|z’; t) to be deterministic or stochas- 
tic, and we consider the former here and discuss the latter in 
Ref. [106]. We take the total cq Hamiltonian to be given by 
A(z) = p?/2m + B(q)H™) with A” = 1 (ata + aa’) 
for the oscillator and A”) = G ©’) for the qubit. 

To model a Stern-Gerlach experiment, we consider a 
two-level system with states |{),| |) and a cq-density matrix 
which depends on the classical degrees of freedom. 


For simplicity, we consider a linear potential B(q) = Bq 
and a decomposition of the Hamiltonian, Eq. (18), in terms 
of Lindblad operators L; = /o|t)(t|, £, = Voll) (1 
with htt(z) =—h!\(z) = Bq and the others zero. An 
example which builds toward understanding the field 
theoretic case is the harmonic oscillator with Lindblad 
operators @ and a". This is more interesting in the sense that 
the pure-Lindbladian term raises and lowers the state of the 
quantum system. In the qubit system, this corresponds to 
choosing Lindblad operators Li =o 14 = 
Vo|*t)()| with diagonal couplings W*(z|z— A) = 
SagW%(z|z — A), or even a single Lindblad operator 


a(q, P) 
u\(q. P) 


uy(q, P) 


law) se 


P(q, P) = ( 


be G7 vey, Here, we just present the more trivial first 
case decomposition; the other cases, including nontrivial 
B(q), are studied both analytically and numerically using 
an unraveling approach in Ref. [106]. We take a simple 
realization of the discrete dynamics, that of Eq. (68) for 
small z, 


(a-« 


S(t) (NA(@. p+ eB)It)(4 


ant 
op 


0 


(78) 
I 
0a(q.P)__4. p 0a(q,p) 
=-—2iB fe wees 
- iBqwa(q, p)—— 7 a(q.p). (80) 
Equation (80) has a solution 
P —2iBat(qg-2)-2! 
aap) = (4-22 )e 2iBot(q-Iy)—e, (81) 
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where p(q — pt) is any normalized function of its argument. 
We see that the coupling of the qubit results in decoherence, 
since the off-diagonal elements decay exponentially fast. 

More than that, the wave function collapses into a 
definite state of being in the up or down state. This we 
see through the equation for the diagonal terms. Comparing 
Eq. (79) to Eq. (12), we see that at a rate of w/z the system 
undergoes jumps in momentum by +7B depending on its 
state. Since we can monitor the classical degrees of freedom 
without disturbing the system, measuring the change in 
momentum will uniquely determine the value of the spin. In 
this sense, the interaction leads to more than just 
decoherence, but an objective collapse of the wave function 
[69,100]. While the basis to which the quantum state 
collapses can appear arbitrary in spontaneous collapse 
models, here it is determined by the form of the CQ 
interaction which will be tightly constrained by the require- 
ment of reproducing Einstein gravity in the classical limit. 
More generally, if the backreaction force Vh® is unique 
when diagonalized, then monitoring the classical system 
uniquely determines which Lindblad operators were 
applied to the state. Thus, an initially pure quantum state 
remains pure conditioned on the classical degrees of 
freedom. In the continuous realizations of the theory, this 
has since been shown to be generically true in Ref. [29] 
with Layton and Weller-Davies, provided an inequality we 
call the decoherence-diffusion trade-off is saturated. 

The solution (a Poisson distribution) for the diagonal 
terms of f in this and a more general setting is given in 
Ref. [106], but the evolution can be understood in broad 
terms by inspection. The particle is making a momentum 
jump of +7B at a rate of w/z, so on average it has an 
acceleration of +wB exactly as we expect for a Stern- 
Gerlach device. We can imagine that at each df, there is a 
probability of the particle undergoing a jump in momen- 
tum, so by the law of large numbers, we expect the particle 
to have made on average wt/z jumps, normally distributed 
with a variance of 2[1 —(@/z)]. If we Taylor expand 
Eq. (79) 


du(q.p) __ pOu(q.P) oB ou(q, P) 
ot m Ox Op 
| Bt Ou, P) eee (82) 
2 ape” 


(where uw is uy or u, for +, respectively), then to first order 
we have a particle undergoing acceleration twB depending 
on its spin, as expected, with a diffusion term at higher 
order. However, care must be taken, as if we truncate the 
Taylor expansion, the density matrix can become negative. 
Equation (78) is completely positive however, and can be 
understood as follows: An initial pure state in superposition 
will initially have a constant momentum, but will then 
undergo jumps in momentum of finite size 7B, and at a rate 


PHYS. REV. X 13, 041040 (2023) 


w/t. This causes the state to collapse to the up or down 
state, at a rate w/t. 

If we trace out the qubit and look at the average value of 
the particle’s trajectory, then for a uniform superposition of 
up and down spin, it will have zero average acceleration, 
but if we look at the full state, the particle undergoes 
acceleration (albeit stochastically) in a direction which 
depends on the value of the spin. We thus see that the 
evolution of the particle corresponds to the left panel of 
Fig. 1, not the right one. 


C. Postquantum field theory 


In order to consider applying our master equation to field 
theories, we must first ensure that it can be made Lorentz 
covariant. Although it was initially thought that such 
dynamics are impossible [137], GKSL equations which 
are covariant under the proper orthochronous Lorentz 
group were introduced in Ref. [138]. This was further 
taken up in Refs. [65,66] (see also Ref. [139]). A Lorentz- 
invariant Schwinger-Keldish action for open quantum 
scalar field theory was used to show that the GKSL 
equation for scalar field theories is renormalizable [140]. 
The formalism used there is manifestly Lorentz invariant. 
There are however subtleties that have since been clarified 
with Weller-Davies in Ref. [141]. The full master equation 
of Ref. [140] is not completely positive, but one can find 
Lorentz-invariant master equations which are CP, Lorentz 
invariant, and renormalizable, as well as master equations 
which are CP, renormalizable, and Lorentz covariant, but 
not Lorentz invariant. The Lindbladian corresponding to 
the discrete realization presented in Sec. V is covariant but 
not invariant, while the Lindbladian of the continuous 
realization is invariant but not covariant. This apparent 
paradox is resolved by showing that transformations of 
space-time are not generated by a unitary, but rather, by the 
semigroup inherited from the Lindbladian, and thus, there 
are different notions of invariance and covariance [141]. 
Here, we focus on covariance conditions, and address the 
subtleties in more detail in Refs. [141,142]. 

In order to make the dynamics here Lorentz covariant, 
one can begin by choosing Lindblad operators and cou- 
plings, such that the right-hand side of Eq. (14) transforms 
like 0/dt, so that both sides of the equation transform in the 
same way under a Lorentz boost. This has been referred to 
as a minimal Lorentz-covariance requirement [137]. We 
first construct a field theory version of Eq. (14) and then 
discuss how to make it covariant under the proper ortho- 
chronous Lorentz group. 

To obtain the field-theoretic version of the master 
equation, we take the Lindblad operators L,, to also depend 
on the spatial coordinate x and in particular, take them to be 
local field operators £,(x). We can think of them as 
carrying a double index a, x. This requires the couplings 
W to carry not just an index a@f but also x, y. Then, instead 
of a sum or integral over just a, #, we also have an integral 
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over x,y, and write the couplings as W7*(z|z’; x, y) with the 
phase-space degrees of freedom z being local fields at each 
point x. 

Cluster separability [64] demands that the coupling 
constants W?(z\|z’;x,y) be delta functions in x, y [60], 
and that the coupling to phase-space variables z be to the 
local classical field z(x). A violation of cluster separability 
can result in violations of Lorentz invariance, with corre- 
lations or entanglement being created over spacelike sep- 
arations. However, the creation of correlations need not 
violate causality or allow for superluminal signaling [65]. 
Thus, while it is desirable to have W°’(z|z’; x, y) be delta 
functions, for greater generality we allow the coupling 
constants W°*(z|z’; x,y) to have support at finite |x — y]. 
| 


PHYS. REV. X 13, 041040 (2023) 


Should these couplings not fall off fast enough with 
increasing |x—yl|, the theory is likely to run afoul of 
experiment. For the cq state at z, we continue to write it 
as P(z;t) with the understanding that it lives in Fock space, 
and is an operator-valued distribution functional of the 
classical fields z(x) over the entire manifold of points x. 
If we consider the points x to be on a lattice, with the 
continuum obtained in the limit that the lattice spacing goes 
to zero, then we can think of the measure Dz over phase 
space as II,dq(x)dp(x), and in the continuum limit, the 
master equation is given by a functional integral. 

If we take the Lindblad operators to be local field 
operators £,(x), then the field-theoretic version of 
Eqs. (14) and (17) is 


Xs: ) _ _ift(2),a(z:8)| + / dxd yDAW# (z\z — Asx, y)La(x)p(z — A; 1)E}(y) 
— 5 Wi (esx, EHO) bala), D(z} (83) 
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One can also consider master equations which contain multiple jumps at higher order, 1.e., 


a 


[Pade rdr,...wePrerte (al —A3X1,X2,...) °° Laas: (x,)A(z — A; OL, (x1) £5, (x2) vee, 


but we do not do so here. 

In flat space, one can consider as an example the scalar 
field discussed later in Sec. V, and take the pure-Lindbladian 
coupling term W%*(z; x, y) to be close to a delta function and 
proportional to 1/z. If we wish to reproduce Hamiltonian 
dynamics in the classical limit, then we generally take the X a 
term to transform just as the Poisson bracket would for a 
Lorentz-invariant classical field theory, since the only differ- 
ence between it and {A(z), A(z; t)} is operator ordering. 

As an example, we could decompose the Hamiltonian of 
the field theory in terms of the Lindblad operators used in 
Refs. [66,138,140] 


A(x) = [asters bf)Lals) (85) 


with h?*(z;.x) coupling the quantum field to a classical field 
z(x), which in the case of gravity, we can take to be the 
3-metric and its conjugate 3-momentum. The couplings 


XP (zg3x) = (86) 


then reproduce the Hamiltonian backreaction of the quan- 
tum fields on the classical ones. If one wants to be explicit 
about W°?(z|z’; x, y), one can use the realization presented 
in Sec. V. 

In Ref. [71], we show that depending on the decom- 
position A%?(z;x), there are discrete dynamics or 
continuous dynamics. In particular, it is shown that the 
Kramers-Moyal expansion in Eq. (84) never terminates and 
leads to finite-sized jumps in phase space unless the only 
nonzero D%*(z; x,y) is D2°(x,y), and the only nonzero 
Xz; x) is X2(z; x) and X(z;x), and the higher-order 
terms all vanish. In other words, the continuous master 
equation is obtained by restricting the decomposition of 
Eq. (85) to include only terms 
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A= / dxh*(z;x)1£.4(x), (87) 


while otherwise, one has discrete dynamics. 

The positivity conditions in the field-theoretic case 
are analogous to the finite-dimensional conditions of 
Eq. (50). Recalling that for 2n phase-space degrees of 
freedom per point, D?°(z; x, y) isa 2n x 2n matrix over x, y 
and denoting W,,(z;x,y) as the matrix with elements 
we (z3x,y), and X),(z;x) as the vector with elements 
X°4(z;x), the continuous master equation is completely 
positive, provided that the square matrix 


D??(x.y) 
X},(z3.x)d(x, y) 


5(x, y)Xp(z5y) 


=0, (88) 
Wi,(z3x, y) 
which for invertible W,(z;x,y) is equivalent to the con- 
dition that D?(x, y) be positive semidefinite along with its 
Schur complement 


DY (x, y) — X}(z3x) Wal (zx, y)Xa(zsy) = 0, (89) 


where W;,'(z; x, y) is the generalized inverse of W;,(z; x, y). 
The positivity condition for singular W,, (z; x, y) is a straight- 
forward extension of Eq. (89); see, for example, Ref. [143]. 
We have since shown that the condition given by Eq. (89) 
holds for the zeroth, first, and second moments of any 
consistent classical quantum coupling, and provides an 
experimentally testable signature of classical gravity [50], 

Turning back to Lorentz covariance, one might worry 
that in the master equation of Eq. (84), there is a 
dependence on the choice 0,, 


OP(zst) = L,p(z3 1). (90) 


As discussed in Ref. [138], we can consider any four-vector 
a from the future cone F*, i.e., a- a > 0, and ap > 0, and 
consider the dynamics 0,/(z;t) = £,. Indeed, consider M 
an element of the proper orthochronous Lorentz group. Let 
m(a, M) be a representation of the Poincaré group and let 
the evolution of the state j(z;t) be given by a completely 
positive trace-preserving Markovian map €,. The Lorentz- 
covariance condition as defined in Ref. [138] (see also 
Refs. [66,144]) is then 


n(a, M)E,x(a,M)" = Ey, (91) 


for any b € F*. In other words, if we perform a boost from 
the lab frame, evolve the state, then boost back into the lab, 
one gets the same final state as evolving the initial state in 
the lab for the boost time. The open quantum scalar field 
theory we consider in the next section satisfies this property 
when the background space-time is flat, as can be seen from 
the fact that in that limit, it is equivalent to the Lindbladian 
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considered in Ref. [138]. The relationship between this 
covariance condition and Lorentz invariance is discussed in 
Ref. [141]. 


V. POSTQUANTUM GENERAL RELATIVITY 


We now have a procedure to take a Hamiltonian, which 
couples classical and quantum degrees of freedom, and use 
it to generate consistent hybrid dynamics. Let us now do 
that for general relativity coupled to quantum matter via the 
energy-momentum tensor. We do so in the ADM formu- 
lation of general relativity [118,145], whose classical 
formulation we now recall. One starts with a slicing up 
of space-time into spacelike hypersurfaces labeled by the 
parameter 7, and with coordinates x. The hypersurfaces are 
endowed with a 3-metric g,,, with determinant g at each 
point x. We often drop the explicit x dependence from field 
degrees of freedom, since all quantities should be assumed 
to be locally defined, unless they are a functional, in which 
case, this is indicated with square brackets. We put in the 
explicit x dependence when we have two integrals over 
space, or if there could be any ambiguity. The components 
of the 3-metric are the canonical coordinates, and x? their 
conjugate momenta 


qe = JVa(K” _ g’K) (92) 


with K@? the extrinsic curvature and K its trace. 

In general relativity, one is free to choose the timelike 
direction by arbitrarily specifying a lapse N and shift vector 
N*“ at each point x. The evolution of the hypersurfaces is 
then governed by the classical ADM Hamiltonian 


Hapm N,N] i / dx(NH) + NeH{) 


4 ¢ ds X F{(boundary) | (93) 


and we now define the various terms appearing in it. 1°) 


is called the total super-Hamiltonian and is typically written 
in terms of a pure gravity part, and a term which contains 
the matter fields 

Hood) — His) +H), (94) 
The first term is the purely gravitational part of the super- 
Hamiltonian given by 


He) i= Caen ~ g?R (95) 
with R the intrinsic curvature on the 3-manifold analogous 


to a potential term, while the first term is analogous to a 
kinetic term, with G,,., the deWitt metric, 


1 
Gobeil = 50 GacIod + GadIbe — Jabed)- (96) 
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H”) is the matter Hamiltonian density, and we take it to 
also include the cosmological constant term \/9gapAce- 


The total supermomentum H\" is given by 


He He ea” (97) 


with H{2) := —2gacV,x°° the purely gravitational contribu- 
tion, and V, the covariant derivative with respect to the 
Levi-Civita connection for the 3-metric. #1" is the matter 
contribution. The final term in Eq. (93) is the boundary term, 
and we ignore it in the present paper, although the formalism 
can be extended to apply to its inclusion. Finally, it proves 


convenient to define the total gravitational Hamiltonian 
H'®)(N,N] and the total matter Hamiltonian H” [N,N] 
as functionals so that Hapy[N,N] = H(N,N] + 
H(™)(N, N] and 


H gr) [N,N] = [axa lite NeH 2), 


HIN, N] = if dx(NH™ + Ne"). (98) 


With the Hamiltonian now defined, let us turn to the 
equations of motion. The dynamics is given by Hamilton’s 
equations 


one? . \] 
’ = ’ 99 
ot Gab ( ) 


O9ab _ OH apm [N, N] 
ot bn” 


where 6 indicates the Fréchet derivative. In order for these 
dynamics to be invariant under the choice of the lapse N 
and shift N“, we need to impose primary constraints on 
their conjugate momenta 

Py, &0, (100) 
where ~ denotes that the constraints are weakly zero in the 
sense that they vanish only on part of the phase space (the 


constraint surface). To ensure that these primary constraints 
are conserved in time, i.e., that 


Py = 0, Py, = 0, (101) 
we require secondary constraints 
HO x0, HY x0. (102) 


The constraints, once initially satisfied, are preserved by the 
dynamical equations of motion. 

One can write a Liouville equation for classical general 
relativity 
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dp(g, #3 t) 


Ot = {Hapm N,N], p} 


= {HN §],p}, + (HIN. §], ohn 


+ {H™[N,N], p},, (103) 


where we divided up the Poisson bracket into one with 
respect to the pure gravitational degrees of freedom 


{H™ [N,N], P}, 
5H \™) \7 5 
= [a( WN] P 


OJab én - 


6H [N,N] dp 
én? OGab 


i: (104) 


and one with respect to the matter degrees of freedom, 
which for a scalar field with conjugate variables ¢,z4 at 
each point x would be 


{H™ N,N], 0} m 


6H\™IN,N| 6p SH (N,N) 6 
= fas [N,N] 5p _ [ 15°) 495) 
op Ong Ong op 


The Poisson bracket with respect to the lapse and shift and 
their conjugate momenta is zero when the constraints are 
satisfied and is usually not included in the definition of 
the Poisson bracket. When we discuss the classical quan- 
tum theory, it may be useful to also include them in the 
definition of the Poisson bracket. 

When only matter is quantum, the matter Poisson 
bracket of Eq. (105) just becomes the commutator 
~i[A” [N, N],/] which is unproblematic, and just repre- 
sents the evolution of the quantum fields on a spatial 
manifold endowed with an evolving classical metric. This 
is just quantum field theory in curved space. However, the 
Poisson bracket of Eq. (104) is the backreaction of the 
quantum fields on space-time, and this is the dynamics we 
now seek to describe in terms of CQ dynamics. 

We first replace T“” with the stress tensor of quantum 


field theory, A") = ar. We can now write 


HUD — pyle 4 Hm, (106) 
as long as we recall that this object is hybrid and cannot be 
used to derive classical or quantum evolution in the 


traditional way. For the discrete master equation, we 


expand A) in terms of local field operators £,(x), so 


that dropping the explicit x dependence 

HO) = He)1 + nF(g, a) EEL, (107) 
with £)(x) = 1 multiplying the purely classical part of the 
super-Hamiltonian. For convenience, we explicitly separate 
out the pure gravity contribution H'"), and one can then 
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take h° = 0. Although in our derivation of the master 
equation, we take the L,, to be a basis and Le to be traceless, 
one can show as with the Lindblad equation, that any set of 
Lindblad operators can be put into this standard form [71], 
and so we henceforth do not assume that the other L, are 
traceless. We often consider matter which is minimally 
coupled, meaning A does not depend on z‘/ so that we 
can write h°’(g). For the continuous master equation, which 
we consider elsewhere [146], we use Lindblad operators 
Lj, = 6H [N, N]/6g;;, H™ [N,N], and Lo = 1. 

The supermomentum can be treated similarly, where the 
stress-energy term in Eq. (97) is replaced by its quantum 
counterpart 
= A” 


HW” (108) 


leading to 


FO? = 1 + pP@ LiL, (109) 
for the discrete master equation. Although the matrix with 
elements N“ pie is not generally a positive-semidefinite 
matrix, the combination Nh? + N4@ pv is. 

To make things concrete, it is useful to have in mind a 
scalar field, especially since the general case can be treated 
similarly and Lindbladian evolution of ¢* theory has been 
shown to be renormalizable in this case [140], albeit with 
the caveats mentioned in Sec. IV C. We have [147] 


A 1 
Hm = 5 (Gag + Vag" V dV ob i Vm? p(x)? 
~ 4/9Mcc)s 


A” = TV ab- 


(110) 


(111) 


For the discrete evolution, we can now decompose the 
Hamiltonian as in Eq. (107), with the metric degrees of 
freedom included in h**(g, z)., 


h™ := ge, hee = V9; 
Ao om Jog” hd = V9 


with the rest zero, leaving the local Lindblad operators as 


(112) 


A A m 
—_—— Le ee, 
var o= 7? 
. 1 , 
La = 0.0. Ly =v 2A|A) (Al. (113) 


Note that the h°(g, z) do not include the N prefactor, since 
we use it to apply the constraint equation; however, when 
we look at the dynamical equations of motion, we need to 
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include it. We choose a simple form for Ly, but a priori, 
other choices are possible, for example, the shift operator. 
Their particular dependence is not important in terms of the 
broad outline of this discussion, rather, what is important is 


that we can decompose #1") in terms of local field 
operators. The expansion of the supermomentum in terms 
of local operators could be treated similarly, with 


L,=tp ILy=00, Ej)=1, (114) 
Note that He”) has no metric dependence and thus does not 
contribute to the backreaction. 

Although we cannot generally define Fock states locally 
on a Cauchy slice, if we know, for example, that the 
dynamics will lead to a space-time which asymptotically 
has a timelike Killing vector, then we can choose a vacuum 


state and set of creation and annihilation operators 
a 5(9)-45(9), which depend on the metric to expand the 
matter Hamiltonian as 


HON.) = | SFeeonl ap oapla) + a9(0)44 (0), 


where we drop the cosmological constant term for simplicity. 
In a fixed background, for example, the Lindblad equation 
using the Lindblad operators of Eq. (113) are equivalent to 
those of Eq. (115). In this way, one sees the covariance of the 
master equation; however, for general dynamical space- 
times, an expansion of local operators in terms of creation 
and annihilation operators is not possible, and we use local 
Lindblad operators such as those of Eq. (113). 

The dynamical equations of motion are generated by the 
CQ-ADM Hamiltonian 


AapmlN. N= / dx(NA + NH”) (116) 


with A) and HW given by Eqs. (107) and (109). If we 
wish the dynamics to be invariant under the choice of lapse 
and shift, this would place strong constraints on it. It seems 
to require that the stochastic part of the master equation be 
generated by NW%*(z|z— A;x,y), ie., that the model be 
linear in the sense that W*/ (z|z — A;.x, y) can be replaced 
by NW? (z\z—A;x,y) and Wi pe(zlz —A;x,y) by 
ap . : : . 
NaWio(zlz- Asx,y). This is analogous to replacing 
{NH , p} by N{H 9, p} in the classical case, since 
they are equal when the constraint HY) ~ 0 is satisfied. In 
what follows, we take the lapse and shift outside the 
functional derivatives under the assumption that they either 
have no explicit metric dependence or the constraints are 


satisfied. The constraints in this theory will not be those of 
GR, but must be modified. As noted in Ref. [156], it is not 
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strictly necessary that the dynamics be invariant under the 
choice of N and N,, and there are sensible theories which 
break this, such as shape dynamics [157-159] or unim- 
odular gravity [160-166], but it is desirable to have some 
local symmetry, even if it is not full diffeomorphism 
invariance. 
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We do not restrict ourselves to a particular realization of 
w*? (zlz—A;x, y) and W#?(z|z — A; x, y) at this point, but 
we assume that the realization satisfies this linear condition 
and then present realizations which have this property. We 
thus have the master equation 


PGi) _ [axe t) A(z: De fee nee (™ IN, N], A(z; 1)] 
ais [ exaypa [ww Wil (zlz— Asx, y)La(x)p(z — A)L} 
if [ exavDa now “h(zlz— Asx, y)Lq(x)A(z - 


where we write z for the phase-space degrees of freedom 
when g,»),2°" are too cumbersome and have absorbed the 
bare term 4° (z;x,y) into W*"(z|z—A;x,y). The fune- 
tional integral DA is over the phase space. Here we take the 
lapse N and shift N“ outside the Poisson bracket, under 
the assumption that the constraints (suitably modified from 
the deterministic theory) are satisfied [156]. 

The first two terms of Eq. (117) correspond to the 
deterministic classical evolution associated with the 
Lindblad operators Ly = 1. The dynamics associated with 
this classical part could be stochastic, and could a priori 
include jumping or diffusion terms as well as friction terms. 
We take such terms to be included in W°(z|z — A;x). 
Because they are associated with Lo, they do not need to be 
positive. 

The next term is the commutator of the state with the 
matter Hamiltonian describing how quantum matter 
evolves in a background space-time. Finally, we have 

| 


+ [ axaywDa| wy Wir (zlz — Asx, y)La(x)A(z — As t)L}(y) = 5 Wi (Eh .Lu (x), A(zst)}4 


This should be thought of as the fundamental dynami- 
cal equation, but we can perform a Kramers-Moyal 
expansion on it. Since we expect all moments to appear, 
the moment expansion will include terms such as 
(5h” /5gav) [6p(z; t)/Sx’| since it transforms like a den- 
sity. This moment at zeroth order in 7 is needed to recover 
the classical limit of general relativity. We now require a 
positive transition matrix W*"(z|z— A;x,y) [except for 
the W9(z|z — A;x, y) component], and that each moment 
transforms as a density to ensure invariance of the 


= | as|n (rte, pleat} + NCH plea} - 


A)L ly 9) 5 Na Wee (23x )(EfOILaC) AD}. (117) 


the interaction terms between gravity and matter which 
must be stochastic. Two particular choices of W°(z; x, y) 
are given as Eqs. (120) and (125), and in both cases, the 
integral over A can be removed. As shorthand, we denote 
Eq. (117) by 


Op(g, mt) 


AE Aaow PCG. 5) 


= ', dx[NLp(9, 231) +N°L,p(9,231)|. (118) 


Because H! ”) has no metric dependence, the evolution along 
N“ can be pure Lindbladian or purely deterministic. In 
Ref. [156], we argue that diffeomorphism invariance requires 


the latter, i.e., that wai (z|z —A;x,y) = Wales y) = 0, 
Equation (117) then simplifies to 


if [N,N], A(z) 


(119) 


I 
dynamics under spatial diffeomorphism. An example of 


such a realization for the discrete master equation is given 
in Ref. [156]; namely, for a, 6 # 0, 


/ dyDANW# (z|z—A;x,y)L_(x)a(z— Ast) Li(y) 


a 


="exp tar) era) ex) -} 


x ht (z;x)Lg(x)p(z; t) L(x x), 


a 


(120) 
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where exp|M]¥ is the components of the matrix exponen- 
tial of the matrix M and it is purely local in x. Here, 
M has components Mf = [c(h~'),;(z){A/*(z),-}] and 
W? (z;x,y) =(1/7)h%(z;x)5(x,y), and gives a Kramers- 
Moyal expansion up to constants of 


/ Wi (cl — A)p(<— Ast)DA 


dh (z) Sp(zst) dh? (z) 6p(z; 1) 


1 7 
= a (Gav)P(z; t) 4 


OGab bx? bn” OGab 
& p(zst) 
ae tDrbed.ap a is ee 2 121 
(9) bn bn ( ) 


where the diffusion terms D4“’°-74(g) and all higher-order 
moments are tensor densities of the appropriate weight to 
contract with the derivatives acting on the density matrix. 
The absence of functional derivatives with respect to g 
higher than 6f/6g,, means that this is a Brownian 
realization. As a result, Eq. (92) is unmodified, and we 

| 


Sh”? . Sp(z3t 
+ [aw Be ae 
O9ab on" 


Let us now see that in the classical limit, we recover 
general relativity to leading order. For convenience of 
presentation, we consider the case where the quantum 
system is uncorrelated with the classical system. The more 
general case is akin to Eq. (23). Recalling Eq. (9), that 
p(z;t) = trp(z;t), we perform the trace on the quantum 
system in Eq. (123). Dropping terms of second order and 
higher we are left with 


op(z; t) 
at 


+ / dx/gttN {Hp} + (124) 


which, when the constraints of GR are satisfied, is the 
classical Liouville equation for the gravitational degrees of 
freedom of general relativity sourced by quantum matter. 
The third term gives the backreaction we expect. We thus 
see that in the classical limit we recover Einstein’s 
equations. To leading order, the evolution of the quantum 
fields is seen from Eq. (123) to correspond to the ordinary 
commutator terms (quantum field theory in curved space) 
plus a Lindbladian term. 


ep. ore ) 
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have only diffusion in 2”. Higher-order functional deriv- 
atives in 2”? can be suitably smeared to control divergen- 
ces. If we want to ensure that our evolution equation is 
invariant under canonical transformations, we can express 
it in terms of Poisson brackets, with the second-order term 
being written, for example, as {J,, {J;, A(z; t)}}. It may 
also be necessary to include the lapse, shift, and their 
conjugate momenta in the Poisson bracket, which is a 
realization we pursued in Ref. [146]. 

There is more freedom when it comes to W9(z|z— A 
x,y) which could include dispersion terms, non-Brownian 
terms, or even friction terms such as those of the form 


by” p a 


on “sab ? (122) 


with 7° a local function of the metric. However, gauge 
invariance places strong restrictions on these terms, and 
here we leave them unspecified and denote them by the 
differential operator T°. The resulting evolution law is 


=f ax | WH ple 1} + NaH, a(zst)} — LH IN, Nal: al 
+ faxtnnit(e) [Lapteie) —F(EpLaatesn)).| + f axnr pie 


& 


<——sRoaa 123 
bn? qed ( ) 


A realization which gives the continuous master equa- 
tion will be discussed in Ref. [50] with Sparaciari, Soda and 
Weller-Davies and Ref. [146], and is given by 


op Z = 
am f aan tneay + f dxne(ta.p} AN. A 
1 A 
45 f aan {He.9)} — (9.70) 
1 
fe 5 | aedyNiaec gixX,y ym pire] 
1 Q A 
45 f axdyn {DP x.»)Jalgex) Adp(4:9)-D}} (125) 
when the constraints are satisfied. Here, £/ := 


(sH™/ 69i;), L := H, and the pure-Lindbladian cou- 
plings 2,,¢c(g;x.y) consist of submatrices with elements 
Aijut = Do,ijnt(Gs%.¥)> Aijnn(Gs%.¥)> Annki(Gs%sY), and the 
constant Annnn(g;x,y). Here, greek indices are akin to 
space-time indices y =n, i, and we use n instead of 0 
since 0 already corresponds to the 1 Lindblad operator. 
The condition for complete positivity can be found from 
Eq. (89). Denoting the generalized inverse of the matrix 
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Do(g;x, y) as D5!, we require that the matrix with elements 
ijkl ap 
Dy (x.y) and — D3" [6T(93x)/89:j) 8 (95 ¥)/69ui] — 
(Do')'"(x,y) be positive semidefinite as a matrix and 
kernel to guarantee complete positivity of the master equa- 
tion. D* (x,y), Ja(N.N, 93%); and Dojjni(N.N, 93%, y) 
are chosen so that the realization is linear in the lapse 
and shift. If we take the J, to be local in the metric, 
the simplest case is then J = \/g, so that the diffusion 


term is 5 f dxdyN(x)Do(x.y) /9(a) V9) 9" (a) 9" (Y)[50/ 
dn") (x)d2*'(y)] +--+, but there are other more complicated 
realizations which also recover the dynamical Einstein 
equation as a limited case (by adding friction and additional 
diffusion terms, for example). The additional Lindbladian 


couplings corresponding to the Lindblad operator 7") are 
not needed in order to reproduce general relativity in the 
classical limit or to preserve positivity, but instead appear to 
be required by gauge invariance. The constraint equations for 
this realization are derived in Ref. [146] by following the 
procedure outlined in Ref. [156]. Alternatively, one can 
include the lapse, shift, and their conjugate momenta as 
canonical degrees of freedom and thus part of the Poisson 
bracket in Eq. (125). 

This realization has the clear advantage of being con- 
tinuous, while the former realization has the advantage of 
being linear in the Hamiltonian. In terms of the evolution of 
the quantum system, for each classical space-time mani- 
fold, the quantum system undergoes both decoherence, as 
well as unitary evolution generated by the Hamiltonian, and 
quantum jumps which are accompanied by a backreaction 
on the space-time. The diffusion terms imply that the 
solution of Eqs. (123) and (125) [or alternatively, Eq. (125)] 
should be a distribution which has a variance about the 
solution to Einstein’s equations when we perform the trace 
over the quantum system. This presents a testable deviation 
from classical general relativity, potentially observable 
when the other terms are small, such as at low acceleration 
[167]. The absence of such terms is likely to falsify this 
model. On the other hand, since diffusion of the metric can 
result in stronger gravitational fields when we might 
otherwise expect none to be present, it raises the possibility 
that diffusion may explain galaxy rotation curves [169] and 
galaxy formation without the need for dark matter. 

In addition to the dynamical equation of motion, 
Eq. (117) or Eq. (123), we need to impose the constraints 
of general relativity. Let us recall how this is done 
classically. There, the Hamiltonian and momentum con- 
straint of Eq. (102) follows from demanding that the 
Lagrangian be invariant under arbitrary time reparametri- 
zations and spatial diffeomorphisms. This leads to the 
primary constraints on their conjugate momenta, py = 0 
and py = 0. Preservation of the primary constraints is 
implemented in the ADM formalism by requiring that 
Py =O and py« =0 (the supermomentum and super- 
Hamiltonian constraints). These are just constraints on 
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the classical phase space, and so we can constrain p to 
lie on the surface in phase space where the constraints 
are zero. 

Were we to impose the constraint in an equivalent 
manner here, it would amount to demanding that El csi 
be invariant under the action of Py and Py . This would 
lead to a Hamiltonian constraint of the form 


He) + 162GH™ 20 (126) 
at every x which would be impossible to satisfy, since it 
contains both a classical term 12") and a local quantum 


operator An. However, such a constraint is not needed; 
rather, we should instead ask, what restrictions are required 
such that the equations of motion satisfy the required 
symmetry? In the classical case, we demand that physical 
degrees of freedom do not depend on the choice of lapse N 
and shift N“%, ie., that these are gauge choices. This 
requirement necessarily leads to the constraints of general 
relativity. Likewise here, we should not first demand that 
the constraints be satisfied, but rather, require that the 
equations of motion be invariant under the choice of lapse 
and shift. This turns out to lead to modified constraints 
which are not of the form (126). We discuss how the 
constraints can be generated from gauge invariance of 
the equations of motion in Ref. [156]. Invariance of the 
equations of motion under spatial diffeomorphisms is 
already guaranteed, as the evolution equation involves 
only integrals over scalar densities. 

In order to find a solution to the dynamical and constraint 
equations, we need the equations to be compatible. 
Attempts to quantize gravity have floundered on this point, 
as the constraints are operators and typically do not 
commute, so no simultaneous solution of the constraints 
exist. This is the case in loop quantum gravity, although it is 
hoped that one can find an operator ordering such that the 
constraints close. Likewise, string theory is background 
dependent, which presumably means that in the full theory, 
gauge invariance is broken. While it is hoped that frame- 
works for string field theory could result in a background- 
independent approach, it is not clear whether this can be 
maintained. Here, we also require constraint compatibility, 
but it takes on a different form. Because we have attempted 
to implement gauge invariance on the level of the equations 
of motion rather than as operator identities, the relevant 
notion of constraint compatibility is not simply the 
commutation of operators, but whether the evolution 
operators and constraints and generators “commute” with 
one another when applied to the state. In other words, we 
require for example 

LalLOAlzt) — L(y) Lalx)p(zst) 0, (127) 
and similarly, we require that the “commutation relations” 
for all other generators vanish on the constraint surface; the 
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right-hand side must be a linear combination of generators 
and constraint equations. This notion of constraint compat- 
ibility replaces the closure of the Dirac algebra in the 
classical theory. 

This program is initiated in Ref. [156] in the Heisenberg 
representation to find gauge-invariant observables A. This 
constraint structure, while involved, is far simpler than the 
constraint structure of loop quantum gravity or canonical 
quantization, and provides a model in which to understand 
some of the conceptual challenges these theories face. 
While the constraints are more complicated because they 
are implemented on the level of equations of motion, the 
classical nature of the gravitational field means they do not 
suffer from operator ordering issues and are more tractable. 
Nonetheless, it may be impossible to find a realization of 
the theory and a symmetry principle, such that the set of 
constraint equations close. In that case, one could still 
impose the standard constraints of general relativity as 
initial conditions, as is done in the ADM formulation, but 
they may not be conserved in time. One could instead 
impose the standard constraints as expectation values, but 
again, it is possible these may also not be conserved in time. 
This would likely mean the theory was not invariant under 
diffeomorphisms of the (3+ 1)-dimensional space-time, 
but it would not render the theory inconsistent, since we 
formulated the theory here, not as a 3 + 1 theory of space- 
time, but rather, as a Hamiltonian theory of the spatial 
metric. The evolution of matter is governed by the 
Hamiltonian and Lindblad operators, which are functionals 
only of the 3-metric. To be consistent, the theory needs to 
be invariant under spatial diffeomorphisms, which it 
manifestly is. However, although it would be a consistent 
theory, it would be unclear how to justify imposing the 
constraint as a mere initial condition. So while the theory 
may be consistent, it would require further motivation. 

An alternative approach can be pursued now that 
classical quantum path integral methods have recently 
been developed [170]. Since the path integral is computed 
in configuration space, this allows for the construction of 
manifestly diffeomorphism-invariant theories. These 
have since been proposed together with Weller-Davies 
[142]. One, based on the trace of Einstein’s equations, is 
different from Einstein’s general relativity, but does at 
least demonstrate that diffeomorphism-invariant classical 
quantum theories are possible in principle. A theory 
which contains all components of the Einstein equations 
is also proposed, but it is not yet known whether 
dynamics which deviate from general relativity are 
suitably suppressed. 


VI. DISCUSSION 


The general form of hybrid classical quantum dynamics 
is presented here as Eq. (14), which includes discrete 
evolution and a continuous form [71]. Based on this, we 
exhibite a theory of classical general relativity coupled to 
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quantum field theory which is invariant under spatial 
diffeomorphisms. The theory is consistent: The dynamics 
is completely positive, norm preserving, and linear in the 
density matrix. We specify the dynamics up to a choice of 
the positive-semidefinite matrix W?"(z|z— A;x,y). The 
choice is constrained by gauge invariance, and thus must 
satisfy further criteria, as discussed in Ref. [156]. Two 
particular choices of W%"(z|z — A;.x, y), are the model of 
Eq. (120) and that corresponding to Eq. (125). These are 
further considered in Ref. [156] and Refs. [50,146], 
respectively, along with the constraint equations which 
correspond to them. In Ref. [71], it is shown that the first 
realization leads to discrete jumps in phase space, while the 
second realization follows from demanding that the dynam- 
ics be continuous in phase space. The first moment in their 
Kramers-Moyal expansion corresponds to general relativity 
in the ADM formulation, and thus both recover classical 
general relativity as a limiting case. A proposal based on the 
path integral is introduced in Ref. [142]. 

We require that the metric remains classical even when 
backreacted upon by quantum fields, and this has two 
important consequences. First, the interaction is necessarily 
stochastic, which is important, since a principle motivation 
behind the theory are attempts to resolve the black-hole 
information problem. That there exists a fully quantum 
theory of gravity has been an argument against theories 
which have fundamental information loss, in the sense that 
given any such dynamics one can always find a purification 
of the state onto some environment, such that the full theory 
is unitary. In this case, Coleman has shown that when we 
trace out the environment, the dynamics is equivalent to a 
unitary theory with unknown coupling constants [62]. This 
renders the information loss trivial: In a universe with a 
finite number of random coupling constants, an observer 
has only to perform a sufficient number of experiments to 
determine these coupling constants, after which time the 
theory will behave arbitrarily close to a unitary one from 
the point of view of this experienced experimentalist. 
However, the argument in Ref. [62] need not apply here. 
Any purification which keeps some degrees of freedom 
classical would need to obey an unnatural condition, 
namely, that it remain of the form 


Pas = >_ 


with system A encoding the classical degrees of freedom 
and system B being the purifying system [171]. For these 
degrees of freedom to remain classical, one further needs to 
impose the condition that one can measure only in the basis 
where |g, z), is diagonal. If one lifts this restriction on 
measurements, or allows dynamics more general than the 
one preserving the form of the state of Eq. (128), then it is 
far from clear that the dynamics would still be completely 
positive and trace preserving. The dynamics here needs 
only to give positive probabilities on measurements which 


g; 1) 4 ® 


gy; 1) ® [Pmatter (9; T)) AB (128) 
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commute with the basis representing the classical degrees 
of freedom. It does not need to give positive probabilities 
for all measurements. Thus, we cannot conclude that this 
theory has consistent dynamics on the purified state. Or to 
put it another way, it is unlikely to be the case that the 
dynamics considered here can be viewed as unitary 
quantum dynamics with a local environment which is 
traced out. The fact that the quantum system couples to 
a classical one, thus provides a potential loophole to the no- 
go theorem of Ref. [62] ruling out nontrivial information 
destruction, while as argued in Sec. I, information destruc- 
tion provides the loophole out of the no-go argument of 
Feynman and Aharonoy against quantum systems interact- 
ing with classical ones. The two no-go theorems annihilate 
each other. 

There may appear to be some tension between the claim 
that the theory may evade the no-go argument of Coleman, 
and the idea that the theory presented here could be an 
effective theory which emerges as the classical limit of a 
fully quantum theory. However, this is not necessarily in 
contradiction. In Ref. [176], we consider fully quantum 
descriptions, which effectively act like a classical quantum 
description. The classical system corresponds to a quantum 
system decohered into a family of coherent states, and for 
some potentials, this description is valid for long times. 
Since this description is fully quantum, Coleman’s argu- 
ment would imply that it corresponds to a unitary theory on 
an enlarged Hilbert space with unknown coupling con- 
stants, and that there exist measurements which would 
allow an experimenter to determine these constants. 
However, in the classical quantum description of this fully 
quantum theory, we do not allow arbitrary measurements 
on the “classical” system—only measurements of g and p 
are allowed. Coleman’s argument allows the experimenter 
to perform measurements in an arbitrary basis, and it 
appears unlikely one could use fixed measurements to 
determine the coupling constants. 

The second consequence of having a quantum classical 
interaction is that the measurement postulate of quantum 
theory is no longer needed. One has fundamental 
decoherence of the quantum field by the classical metric. 
In this view, classicality is not an emergent or effective 
property of large quantum systems, but rather, fundamental. 
In the discrete theory discussed here, h®% is diagonal and 
unique for each Lindblad operator. As with the Stern- 
Gerlach example, the accompanying jump in phase space 
unambiguously determines which Lindblad operator was 
applied to the quantum state. If we know the initial pure 
state of the quantum system, then by monitoring the 
classical system, we know which sequence of Lindblad 
operators were applied to the quantum state and at what 
times. As such, the quantum state conditioned on the 
classical degrees of freedom remains pure. Likewise, in 
the continuous case, we have since found that there is a very 
natural class of theories in which the quantum state remains 


PHYS. REV. X 13, 041040 (2023) 


pure conditioned on the classical trajectory [29]. This is 
ironic given the motivation of obtaining dynamics which 
allowed for the destruction of information in black-hole 
evaporation. The dynamics here, while stochastic, can leave 
the quantum state pure. It is the classical degrees of 
freedom which gain entropy. This suggests that in black- 
hole evaporation, the associated entropy is the entropy of 
the classical space-time. Since the formalism introduced 
here allows one to study the backreaction of matter fields 
on the gravitational field, there is hope that we can better 
understand black holes, as well as other systems where 
gravity and quantum effects are important. 

Let us next address a number of open questions and 
challenges that the theory faces, starting with the objection 
to theories of information destruction due to Banks, Peskin, 
and Susskind (BPS) [60]. BPS argued that a local inter- 
action which leads to decoherence will lead to anomolous 
heating to such a degree that it will contradict observation. 
Although the BPS argument was not made in the context of 
classical quantum theories, it does apply to Lindbladian 
evolution of the quantum state 6 with Lindblad operators 
being local field operators, and our theory does contain such 
terms. In particular, BPS considered an evolution equation 


(0) _ pt 6(0) 
* at dxdx' D6? (x — x(t, (x)6(1)L5(x’) 
(129) 


and couplings D@’(x — x’). The Hamiltonian of the matter 
fields A”) will not be conserved, and when D6? (x — x’) has 
short range (decays rapidly with |x — x’|) there may be 
significant anomalous heating. In essence, local operators 
L,(x) do not commute with the Hamiltonian, and in the 


limit where D¢?(x — x’) is a delta function, the average rate 
of energy production contains terms which diverge. This is 
easily seen in the Heisenberg representation, since when 


DP (x — x’) is symmetric in x, x’, and a, f, 


Ey(m) 
=. ii dxdx’ Dy(x—x')[4(x), (A, 24x]. (130) 
If the £,(x) contain only a finite number of derivatives of 
local field operators, then the commutator gives terms 
proportional to 6(x — x’), or kth derivatives V*5(x — x’). 
These terms will diverge if Do(e) diverges as e > 0. For 
example, if Do(x — x’) « d(x — x’). 

This anomalous heating should trouble us if it contra- 
dicts observation, but it need not trouble us on a funda- 
mental level. If the dynamics is Lindbladian, there is no 
reason to expect energy conservation; the dynamics is not 
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unitary and so Noether’s theorem does not apply. The 
generator of time translations is not the Hamiltonian but the 
Lindbladian [177]. In fact, there is no way to canonically 
define energy, since by choosing a different basis of 
Lindblad operators, the Hamiltonian will change. In the 
case of classical quantum systems, the notion of energy 
becomes even harder to define, since the Hamiltonian is a 
hybrid object. On top of this, in general relativity there is no 
local definition of energy or conservation law for energy 
integrated over an infinitesimal spatial surface [180]. From 
a physical point of view, this is arguably the more relevant 
quantity, and so the question of energy conservation is not a 
matter of principle, but rather, consistency with experiment. 
One must ensure that the vacuum is stable enough that it 
does not contradict experimental observation. 

The vacuum can be made stable by turning down the 
Lindbladian coupling so that Do(e) decreases as € > 0. 
Since dH") /dt goes as Do(e), we can regularize it by 
discretizing space-time, and scaling Do(e) with e. In the 
pure-Lindbladian case, there is no reason not to take the 
Lindbladian coupling to be arbitrarily small, in which case 
the effect is weak in flat space and only strong when the 
3-metric becomes singular, since we expect Do(x — x’) to 


contain terms proportional to \/g(x) and \/g(x’). Seen in 
this light, the problem raised by BPS is merely one of 
having to introduce another scale. 

A similar issue occurs in spontaneous collapse models 
[67-69], but in and of itself, this does not present a 
problem, since again, one can tune down the 
Lindbladian coupling Do(e) as small as one wishes. 
Those who advocate for spontaneous collapse models 
typically require Do(x — x’) to be large enough to account 
for the level of decoherence we see around us. This 
requirement on Dy is not needed here for three reasons. 
First, environmental decoherence is expected to account for 
most of the decoherence we observe. Second, the collapse 
of the wave function is a byproduct which occurs only 
when sufficient correlations has build up between the 
gravitational field and the quantum state together with 
its correlated environment. And third, decreasing Do(e) 
requires more gravitational diffusion [50], which causes 
additional decoherence once the gravitational field is 
integrated out [30,92]. The heating effect of this secondary 
decoherence is a source of concern, but appears to be 
milder than in BPS [30,50], although whether this is within 
experimental limits requires further investigation since at 
short distances, we are operating in a regime where 
relativistic effects become important. 

The question of whether the couplings can be tuned 
down as we decrease ¢ is related to the question of whether 
the theory is regularizable [170]. In the case of a scalar field 
with a single Lindblad operator #(x), the Lindbladian 
dynamics is both renormalizable and Lorentz invariant 
[140,141], and we can choose the coupling constant to 
generate as little anomalous heating as we wish, since we 
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do not require the amount of entropy production to be large. 
Understanding whether the theory can be regularized, 
remains an important open question for the program 
presented here. The question is complicated by the fact 
that relativistic corrections in the gravitational field become 
important at short distances [168] and by the fact that some 
aspects of the heating appear to be gauge artifacts [181]. 
Other strategies have been considered, but they may 
make the theory less compelling. One may be able to make 
the anomalous heating small by introducing a small 
violation of locality as has been attempted in the case of 
Lindbladians in the context of field theories [66], or as is 
done in spontaneous collapse models where the coupling is 
modified at a length scale at the order of the proton radius. 
It is worth emphasizing that if one reduces heating by 
allowing D#’(x — x’) to have finite range, this does not 
necessarily violate locality. If £,(x) are local operators, 
then it is possible for local observables to still evolve 
locally even though pe (x -— x’). The price we pay if 
D@ (x — x’) does not fall off fast enough with |x — x'| is 
that correlations can be created over finite distances 
[60,65]. While this seems to suggest some correlated 
environment should be evoked, it does not lead to super- 
luminal signaling or inconsistencies. It will however 
naively violate Lorentz invariance, since one has spatial 
correlations but no temporal ones. This could be rectified 
by taking pe (x — x’) to not depend on the spatial distance 
|x — x’|, but rather the space-time distance. For example, 
one could take it to be the Green’s function of the Laplace- 
Beltrami operator. This would then suggest a non- 
Markovian evolution reminiscent of an effective theory 
rather than a fundamental one. Another alternative is to 
modify gravity at short distances. Since gravity is tested 
only to within the milometer scale, such modifications have 
often been considered, but would likely necessitate some 
additional structure such as extra dimensions [182]. 
Whether such measures are required is unclear. In the 
context of the theories considered here, the situation is 
significantly different from the pure-Lindbladian evolution 
considered by BPS. First, the action of the Lindblad 
operators on the quantum field is accompanied by a back- 
action on the classical degrees of freedom, and in particular, 
the gravitational degrees of freedom. While the argument of 
BPS suggests there will be an increase in the expectation 
value of the energy density of the matter Hamiltonian 
density A”), this can be accompanied by a decrease in 
the energy of the classical system due to this backreaction. In 
this case, this is the kinetic energy density G; ,)2'/2"' of the 


pure gravity part of the super-Hamiltonian, #1"). In the case 
of a discrete realization of this theory, such as Eq. (120), the 
action of the Lindblad operator is accompanied by a jump in 
n'J(x). Likewise, the diffusion term of the continuous 
realization, Eq. (125) acts on the kinetic term of H'") (x) 
in the Heisenberg representation producing an energy 
ijkl 


density flux of Dj” (x, x)G;j,)(x) which can be negative 
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and state independent, and also needs to be regularized. For 
example, consider 


| Noon a; 
lim, D(x, x") = Do (git + gig — 261g) 
XA => 


(131) 


with Dy some constant which depends on the cutoff 7, and 
ijkl 


B < 1/3 to ensure D;™ is positive semidefinite. If we wish 


dH'®)(N, N] 
dt 


sy oA = 1 ap 
= {H®") [N,N], Hapm(N. N]} +5 i dxdx' Di" (x, x’) 
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Des x’) to include only the final term in Eq. (131) 
proportional to g'/g*, then we can take D, > D>/f and 
p- -o. 

Using the diffusion coefficient of Eq. (131), we can 
compute the apparent flux of gravitational kinetic energy 
density using the continuous realization of Eq. (125) in the 
Heisenberg representation. Ignoring any boundary terms, 
this gives 


& 


___—___ (8) [N,N 
dr'i (x) da*! (x') Ne) 


- : : i oon ee 
= {H®) [N,N], AN, N]} + [awp.va, (g* gi + gg — 28g 9 Gin) 


= {H®)[N, N], HN, N]} + ( 


where d is the number of spatial dimensions. We assume 


that the lapse and shift NV, N depend only on the metric and 
not its conjugate momenta. The first term is just the change 
in energy due to the backreaction of the matter Hamil- 
tonian, and is compensated for by the action of H'®) [N, N] 
on A’) (N,N] which results in an equal and opposite 
energy flux —{H'®)(N,N], A [N,N]} into the matter 
degrees of freedom. It is the second term which is due 
to the diffusion. Taking #< 1 results in an apparent 
negative flux of gravitational kinetic energy. It can at least 
partially compensate for the heating of the matter field 
given in Eq. (130). The question is a subtle one; for 
example, in the weak-field limit, what initially appears to 
be constant anomalous energy production in the gravita- 
tional degrees of freedom turns out to be a gauge artifact 
[181]. The dependence of Eq. (132) on the lapse suggests 
this might be the case more generally. 

Whether this backreaction is enough to restore equilib- 
rium or slow the production of energy in the matter sector 
of the theory is an open question. It may merely be that the 
total Hamiltonian is conserved on average, with heating 
occurring in the matter degrees of freedom, and a negative 
flux of kinetic energy occurring in the gravitational degrees 
of freedom. To what extent this would be observable has 
since been partially taken up in Ref. [50], where it is shown 
that the decoherence (and heating) in the matter degrees of 
freedom can be made small, but only at the expense of a 
large diffusion in the gravitational degrees of freedom. 
However, because gravity is so weak, the effect appears to 
be within current experimental bounds, although they do 
constrain Di(x, x’). The question then becomes whether 
the anomalous heating and diffusion are within experimen- 
tal bounds, while at the same time, respecting various local 
symmetry properties such as diffeomorphism invariance or 
some weaker version. 


oS) +p, [ asn?, 


(132) 


Related to this is the question of whether the 
Hamiltonian and momentum constraint is preserved in 
time. In the deterministic theory, this says that the total 
energy and momentum density is everywhere zero. These 
constraints are however modified in the theory considered 
here. Since conservation of the constraints in general 
relativity enforces diffeomorphism invariance, the issue 
of diffusion of gravitational kinetic energy and anomalous 
heating is linked with whether diffeomorphism invariance 
or a related local symmetry will be respected. Whether the 
theory satisfies some local symmetry then becomes a 
question of satisfying some potentially modified constraint 
equations, and whether one can satisfy them all at the same 
time. This turns on the question of the compatibility of the 
constraint equations, i.e., whether the algebroid of con- 
straints [183] closes. Of course in the CQ theory, the 
constraint is not of the form Eq. (126), and the full 
dynamics are subtle. A fuller understanding of the algeb- 
roid of constraint equations thus remains an open issue, 
with some initial progress reported in Ref. [156] where 
these issues are further explored (cf. Ref. [142]). 

There is some reason to believe 3 + | diffeomorphism 
invariance or something related to it can be respected, given 
the diffeomorphism-invariant realization of Ref. [142] 
based on the trace of Einstein’s equations, and using the 
path integral formulation of CQ dynamics [170]. In 
Ref. [142], we also present a diffeomorphism-invariant 
theory which gives all of Einstein’s equations, although we 
have not proven that the dynamics preserves positivity of 
initial states. The challenge then is to see if the realization is 
consistent with experimental constraints such as those on 
anomalous heating and gravitational diffusion, and is 
ideally also renormalizable in the matter degrees of free- 
dom. While realizations exist which are renormalizable, 
and realizations exist which have small amounts of 
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anomalous heating, the question of whether a realization 
exists which satisfies all these requirements remains open. 

The classical quantum theory has a number of exper- 
imental signatures, some of which have since been quanti- 
fied in Ref. [50]. It predicts a form of gravitational 
decoherence, which already a number of experiments are 
looking for [184-186]. It predicts stochastic, finite-sized 
jumps, or continuous diffusion of the gravitational field 
which might be detectable in the lab with current technol- 
ogy. The theory also predicts diffusion which might be 
observable in astrophysical systems and in cosmology. The 
theory is highly constrained, especially if 7 the parameter in 
Eq. (17) is taken to be of order h, since increasing the 
coherence time results in larger jumps and more diffusion. 
Finally, the locally correlated realizations considered here 
predict a null result to recently proposed experiments to test 
the quantum nature of gravity [187-190]. 

Here, space-time is treated as fundamentally classical, 
but one could instead view this theory as a candidate for an 
effective theory which results from taking the classical limit 
of the gravitational degrees of freedom of a fully quantum 
theory of gravity. There are a number of ways of taking this 
partial classical limit. The extent and the conditions under 
which this gives the equations of motion presented here is 
taken up in Ref. [176]. We find that care should be taken, 
since an effective theory might violate our assumptions of 
complete positivity at short timescales or of being time 
local. Nonetheless we find that there is a regime, and a 
choice of parameters, such that Eq. (14) describes the 
consistent evolution of interacting quantum systems in the 
limit where a subsystem behaves classically. 

Whether the theory presented here is taken to be 
fundamental or not, it provides a consistent way to explore 
the regime where space-time is classical in a manner which 
does not suffer from the pathologies of the semiclassical 
equations. This has been lacking in our investigation into 
how space-time reacts to quantum systems. This gives hope 
that the ability to consistently study backreaction effects 
finds application on large-distance scales, in cosmology, 
and in black-hole evaporation. 
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APPENDIX A: CONDITIONS FOR COMPLETE 
POSITIVITY OF A CQ MAP 


Here we extend a theorem due to Kossakowski [194] for 
a cq map €(t) to be completely positive on a separable 
Hilbert space; see also Ref. [50] for a more complete proof 
using indicator functions. Let us consider all complete sets 
of projectors P,,(z) on the doubled Hilbert space Hs ® Hs 
and classical degrees of freedom z. 

Theorem 1. £ is a generator of continuous-time, 
completely positive, Markovian evolution of A(z; ft) which 
preserves f dztrp(z;t) if and only if for all sets of 
orthogonal projectors {P,.(z)}: 


trP.(z)Lyz @ 1(P,(z')) 20, Wr#s,zz, (Al) 
trP,(z)Lyy @ U(P.(z’)) 20, Vs,z#7, (A2) 
' trl ® Mey H=0. Va, (A3) 


where LA(z;t) = f dz'Ly p(z'st). 
As a consequence of the above, we automatically get 


tf d2'5(z-2')Ps(z)Ly7 @1(P.(z7)) <0 Vs,z — (A4) 


since the probability conservation condition, Eq. (A3) 
ensures that outward flux from any P,(z) contributes 
negatively to the probability rate. The proof of Theorem 
1 follows from writing 


E42 (6t) = 16(z ~ z) 
m +f 
6t0 ot 


(AS) 


where E(t)p(z;t) = f dz'E,,(t)A(z'st), and using the fact 
that the map €,,, is completely positive, probability- 
conserving map if and only if 


trP.(z)Ej2 @ U(P.(z)) 20 V{Ps(z)} {P(Z) 2% 2 
(A6) 


and 
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S tej @ 1(P,(z')) = 1 V{P,(2)}.r.z. (AT) APPENDIX B: HEISENBERG REPRESENTATION 


Equations (14) and (17) tell us how the classical quantum 
state evolves. As with both quantum and classical mechan- 
ics, there is an alternative “Heisenberg” picture in which the 
state does not evolve but the operators do. One can easily 
derive the time evolution of operators in this representation, 
since it follows from demanding that for all f(z; 1), 


It follows that €,) is trace nonincreasing 


Ejs @U(P,(2)) <1. (A8) 


4, | aettdstesn)isle) = 5 f detrAn(zit\Pu(2). (BL) 
Equation (A2) is an additional constraint on the quantum 
system not found in Ref. [194], although if one considers 
the projector to also be on the classical degree of freedom, | where the subscript S indicates the Schrédinger represen- 
then it becomes subsumed by Eq. (Al). From Kraus’s _ tation, and H the Heisenberg representation. Using this and 
theorem, we know that this implies that A°?(z|z’;r) in Eq. (14) for the evolution of f(z) leads to the Heisenberg 
Eq. (34) needs to be positive semidefinite in a, # forall z,z’. evolution equation for the operator A(z; t) over phase space 


A 


dAg(at) _ An(&)  sac aces) 


dt Ot 
ae ree Tr 
+ i. dhW®(z + Alz)LjAn(z + AOL, - 5 WO (2) {Lh pla Auli i. 


_ 0Ay (zt) 
Ot 
— Lj {h®, Ay(zst)}£] +++, (B2) 


+ ifA(z), Au(z; 1)] + W%(z) |LpAn (zs )L 4 — 


where we use cyclicity of operators under the trace and integration by parts. We include the purely classical term in W°°, but 
for deterministic evolution, it would be —{H(z),A,(z;t)}. For a field theory, and in general relativity, we are usually 
interested in the evolution of diffeomorphism-invariant operators such as 


Ay(z) = f arvaA(eix), (B3) 


where A(z; x) is a local scalar operator. In such a case, we have 


ae _ “Al - (2), Anse: 
+ i; dx [aawen( = & 2x, y)Lp(y)An(z +A, t)£,(x) - 5W(zsx,y{EPO)Lu(2). Ane 0}. 
= Anlst) + i[A(z),Ay(z;t)] + | W?(z;x,y) EjoVAnle t)L,(x) - 5 (Ljo)balx). Ane | ae 
— fsx) Wei), Au (zit) hala) +> a 


Note that care must be taken here. The Heisenberg representation allows one to compute the change of the expectation value 
of an observable A with time. However, if the evolution is nonunitary, then it is not necessarily the case that 
(A"(t)) = (A(t))". For example, if there is dispersion, than typically (2) — ()? will increase with time. Thus, knowing 


A 


A(t) does not necessarily allow one to deduce the expectation value other moments of the operator. 
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